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ABSTRACT 



The lattice Coulomb-gauge hamiltonian is derived from 
the transfer matrix of Wilson's Euclidean lattice gauge theory, 
wherein the lattice form of Gauss's law is satisfied identically. 
The restriction to a fundamental modular region (no Gribov 
copies) is implemented in an effective hamiltonian by the 
addition of a "horizon function" G to the lattice Coulomb- 
gauge hamiltonian. Its coefficient yo is a thermodynamic 
parameter that ultimately sets the scale for hadronic mass, and 
which is related to the bare coupling constant go by a "horizon 
condition". This condition determines the low momentum 
behavior of the (ghost) propagator that transmits the 
instantaneous longitudinal color-electric field, and thereby 
provides for a confinement-like feature in leading order in a 
new weak-coupling expansion. 
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1. Introduction 



Quantum chromodynamics, which appears to provide a satisfactory 
framework for understanding the strong interactions, is expected to have 
two striking features: (i) the force between colored particles is of such 
long range that it confines them, and so (ii) the gluons whose exchange is 
responsible for this force are absent from the physical spectrum. These 
two features, though complimentary, appear difficult to accommodate 
simultaneously in a manifestly Lorentz-covariant theory where one might 
naively expect that the gluon propagator at zero momentum should (i) 
possess a strong singularity as a manifestation of the long-range force and 
(ii) be strongly suppressed because there are no massless gluons.* By 
contrast, these two features are compatible in the Coulomb gauge because 
the long-range force can be carried by an instantaneous static color- 
Coulomb field, while the propagator of (three-dimensionally) transverse, 
would-be physical gluons is suppressed. The instantaneous field is required 
by Gauss's law, which is satisfied identically by the color-electric fields in 
the Coulomb-gauge hamiltonian. The restriction to the fundamental 
modular region, a region without Gribov copies [1], makes the Gauss's law 
constraint more severe. Indeed many years ago Gribov [1] proposed this as 
a confining mechanism. 



The ultimate consistency of a Lorentz or Euclidean-lattice covariant gauge is not in 
question. However the Coulomb gauge presents a simpler picture of confinement and 
allows more extensive calculations than the author has so far achieved in a covariant gauge 
despite considerable effort. 
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We derive the lattice Coulomb-gauge hamiltonian by gauge-fixing in 
Wilson's Euclidean lattice gauge theory. This is a quantized, regularized 
and gauge-invariant theory which one is free to gauge fix. After gauge- 
fixing, the lattice Coulomb-gauge hamiltonian is obtained from the transfer 
matrix of the Euclidean lattice theory by integrating out the variables on 
the vertical links in the limit in which the length of these links approaches 
zero. One finds that the lattice form of Gauss's law is satisfied exactly, and 
an instantaneous long-range confining potential appears in the lattice 
Coulomb-gauge hamiltonian. The gauge-fixed theory inherits the non- 
perturbative regularization of Wilson's theory, so ultra-violet and infra-red 
divergences are strictly controlled. Although cancellation of the ultraviolet 
divergences in the lattice Coulomb gauge is not yet established explicitly, 
this should be a technical problem as long as the cancellation occurs in 
Wilson's Euclidean gauge theory. 

In classical Yang-Mills theory. Gauss's law in the presence of quark 
color-charge density pqu reads 

(D-E)a = V-Ea + gofabcAb-Ec = go pqu^, (1.1) 

where A represents the 3 spatial components of the vector potential, and D 
is the 3-dimensional gauge-covariant derivative. (Vectors are always 
3-vectors.) Gauss's law fixes the longitudinal component of the color- 
electric field. Indeed, if we decompose E into its transverse and 
longitudinal parts. 
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Eia(x) = Eitr>a(x) - Vi(pa(x), 



(1.2) 



where V Etr = 0, and cp is a color-Coulomb potential, then Gauss's law 
reads 

- (D- V(p)a = - go fabc Ab.Etr,c + gg pq^a = gQ pa. (i .3) 

Here a color-charge density pa is defined in which only Etr appears instead 
of E. The operator on the left - where one would find -V^ in an Abelian 
theory - is the 3-dimensional Faddeev-Popov operator M(A) = - D(A)-V. 
Gauss's law determines the color-Coulomb potential (p in terms of the 
remaining variables by means of the instantaneous and Faddeev-Popov 
propagator M-l(A), 

(pa(x) = go (M-1 p)a(x) = go J d3y (M-l)a,b(x, y; A) pb(y). (1 .4) 

Despite its non-invariant appearance, the color-Coulomb law holds in every 
gauge and in every Lorentz frame. The color-electric field 

E = Etr - go V M-1 p, (1.5) 

Ea = Etr>a . go V (M-l)ab (. fbcd Ac-Etr.d + p^^b) (1 .6) 

is expressed in terms of the transverse color-electric field and (the spatial 
components of) A, 
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In the Coulomb gauge, one fixes Ai^(x) = Aitr,a(x). This makes 
M(A) = MCAtr) symmetric. It also fixes AO = - Aq by a time- 
independent constraint. For with E = - dA/dt - DA^, one has 
V-E = MAO = . v2(p, so 

AO = M-1(-V2)(p = goM-l(-V2)M-lp. (1.7) 

In the hamiltonian formulation in the Coulomb gauge, the conjugate 
canonical dynamical variables are Atr and -Etr. The expressions that appear 
on the right-hand side of the last four equations are functions of these 
dynamical variables only. In the quantum theory, one sets Etr = id/dA^^, 
and the continuum quantum Coulomb hamiltonian of Christ and Lee [2], is 
given by 

H = 2-l Jd3x(a-lEtGE + B2), (1.8) 

Here E is as above, £ijkBk^ = ViAjtr'a . VjAitr'a + gQfabcAitr,bAjtr,c^ and 
a = detM(Atr). The energy of the color-electric field (1.6) includes the 
instantaneous longitudinal field produced by the A-dependent integral 
operator M-l(Atr). We shall see that the restriction to the fundamental 
modular region makes this a confining field. 

The lattice Coulomb hamiltonian we shall obtain is of canonical form 
and contains color-electric fields that satisfy the lattice form of Gauss's law 
identically. Indeed, in the coordinate system which we shall use, which on 
each spatial link is. 



5 



t^Ax,i^ = 2-l(Ux,i - Ux,it)traceless, 



(1.9) 



equations (1.2) - (1.6) are unchanged in form in going from the continuum 
theory to the lattice. (Here Ux,i is Wilson's link variable that transforms 
under gauge transformation according to Ux,i gx^Ux,igy, where 
y = x+ei. To avoid possible confusion we emphasize that gauge-fixing is 
done after Wilson's gauge-invariant lattice quantization by compact 
variables.) The only change in formula (1.8) for the continuum quantum 
hamiltonian is that in the lattice Coulomb hamiltonian, the link variable 
Ex,i,a is replaced by 



where the Ja^(A) satisfy the differential equations of the Lie algebra of 
SU(n) on each link, namely [JaC(A)a/aAc, Jbd(A)a/aAd] = fabcj^,d(A)a/aAd 
(and of course the magnetic energy is the plaquette action on spatial 
plaquettes, and the weight o includes the Haar measure). The lattice form 
of equations (1.4) - (1.7), and the Heisenberg equations of motion of the 
lattice Coulomb hamiltonian provide a quantized and regularized 
expression for all the classical Yang-Mills equations of motion. Use of the 
coordinates (1.9) instead of the exponential mapping Ux,i = exp(taAx,i^), 
may be essential to obtaining the lattice Coulomb hamiltonian in closed 
form, which otherwise appears to be known only as an infinite power 
series [3]. 

In the quantum gauge theory, in addition to the lattice Coulomb 
hamiltonian, one specifies a configuration space which is free of Gribov 



E 



x,i,a 



= Jab(Ax,itOEx,i,b, 



(1.10) 
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copies. (For an alternative approach, see [4] and [5].) This is the 
fundamental modular region described in sects. 3 and 5. Section 6 contains 
the crucial observation that the restriction to the fundamental modular 
region contradicts the usual perturbative expansion in powers of go. We 
find however that it is consistent with an expansion in powers of go, in 
which the leading contribution to <M-1> is of order go"l, 

< (M-l)a,b(x, y; Atr) > = go"! u(x-y) 5a,b + 0(1). (1.11a) 

This gives an x-dependent mean-field contribution to the inverse Faddeev- 
Popov matrix itself 

(M-l)a,b(x, y; Atr) = go'l u(x-y) 6a,b + 0(1). (1.11b) 

Thus, by (1.5), the color-electric field acquires a longitudinal part that is of 
order go^ instead of go^, 

E = Etr - V 1 d3y u(x - y) p(y) + 0(go) (1.12) 

which contributes to the zero-order hamiltonian. This is a rather literal 
expression of anti- screening. We call u(r), which is the leading-order 
mean-field contribution to goM-1, the static color-Coulomb potential. It is 
independent of go and Atr, ^nd is determined by a consistency condition, 
eq. (6.22) below. 

The organization of this article is as follows. Section 2 is didactic in 
style, although the results are believed to be new. It addresses the problem 
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of how a restriction on the configuration space, such as the restriction to a 
fundamental modular region A, may be implemented for a system with a 
large number of degrees of freedom, such as a lattice theory described by a 
hamiltonian H. The implementation relies on the close connection between 
a quantum mechanical hamiltonian and a classical statistical mechanical 
system with partition function Z = tr[exp(-PH)] = limM^ooTM, where the 
transfer matrix is given by T = exp(-8H). Suppose that the configuration 
space is restricted by a condition of the form G(A) < 0, where G(A) is a 
bulk quantity. In sect 2 it is shown that in the partition function the 
restriction may be implemented by using the transfer matrix 
Teff = exp[ -£ (H + yG)], and integrating over the unrestricted 
configuration space. Here y is a thermodynamic parameter that is 
determined in the infinite- volume limit by the condition 



where the expectation value is calculated in the system with partition 
function Zeff = limM^ooTeff^. Because of the equivalence between a 
statistical mechanical system and a quantum mechanical system, this 
corresponds to using the effective hamiltonian. 



in the unrestricted configuration space, and calculating <G) in the ground 
state of Heff. Thus the familiar equivalence of different statistical 
mechanical ensembles, such the micro-canonical and canonical, implies the 
equivalence of the corresponding quantum-mechanical systems. We call G, 



limv^ooV-l<G> = 0., 



(1.13) 



Heff = H + yG, 



(1.14) 
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the horizon function, and (1.13) the horizon condition. The method is 
illustrated by some elementary examples. In particular, the eigenvalues of 
the N-dimensional laplacian operator restricted to the interior of an 
N-dimensional ellipsoid are calculated in the limit N ^ ©o. The imple- 
mentation of the restriction of the configuration space by (1.13) and (1.14) 
is correct only in the infinite-volume limit, where the number of degrees 
of freedom becomes infinite. On the contrary, in the works of Liischer [6], 
Cutkosky [7], and van Baal [8] and co-workers, which also use the Coulomb 
hamiltonian, the restriction to the fundamental modular region is imposed 
in the approximation in which a finite number of modes are retained. 

For the case of QCD, it will be found that the thermodynamic 
parameter y sets the scale for hadronic mass, and corresponds in 
conventional perturbation theory to the introduction of a renormalization 
mass Aqcd. The horizon condition (1.13) may be solved for go instead of 
Y, which corresponds in conventional perturbation theory to solving the 
renormalization-group equation for go. 

In sect. 3 we give a precise definition of the minimal lattice Coulomb 
gauge, which is used to fix Wilson's Euclidean gauge-invariant lattice 
theory. It is has the defining property that the link variables U are made as 
close to unity as possible, in an optimal way, on all the spatial links of the 
Euclidean lattice. Some elementary properties of this gauge are derived in 
sect. 3. 
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In sect. 4 and Appendix A, the lattice Coulomb hamiltonian HcouL 
eq. (4.16), is derived from the transfer matrix associated with the partition 
function of Wilson's Euclidean lattice gauge theory, 

Z = j dU exp[ - Sw(U) ] . (1.15) 

The variables associated with the vertical links are integrated out in the 
limit in which the length ao of the vertical links approaches zero, so the 
hamiltonian depends only on the variables associated to horizontal links of 
the lattice. The Heisenberg equations of motion of this hamiltonian, and 
the constraints that are satisfied identically, provide a quantized and 
regularized version of all the classical Yang-Mills equations. 

In addition to the hamiltonian, the quantum theory also requires 
specification of the configuration space. In sect. 5, the "horizon function" 
G(Atr) is defined. It describes the fundamental modular region A of QCD 
in the minimal Coulomb gauge on a periodic lattice in the infinite-volume 
limit by the condition G(Atr) < q . (Its expression on a finite lattice is not 
known.) The original derivation of this function in the continuum [9] and 
lattice theory [10] has been improved recently [11], but it still cannot be 
regarded as rigorous. However, it has recently been discovered [12,13], 
following the work of Maggiore and Schaden [14], that G(Atr) vanishes 
identically, G(Atr) = q for all Atr, on a finite lattice with free boundary 
conditions, as is explained in sect. 5 (although not of course on a finite 
periodic lattice). This is just what is wanted, because in the infinite-volume 
limit of the periodic lattice with the horizon condition, (G>/V = 0, the 
distinction between free and periodic boundary conditions should be lost. 
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Thus the modified hamiltonian Heff = Hcoul + YoG, differs from the 
canonical hamiltonian Hcoul by something that is formally zero. 

The theory described by Heff, eq. (1.14), is not a gauge theory unless 
the horizon condition (1.13) is satisfied, which fixes y to the critical value 
Ycr = Y(go)- In sect. 6, it is shown that the horizon condition contradicts 
the usual perturbative expansion. We therefore introduce the new 
expansion in powers of go, eq. (1.11), and derive the self-consistency 
equation (6.22) which is an integral equation that determines the static 
color-Coulomb potential u. 

In sect. 7, the new expansion (1.11) is everywhere substituted into 
Heff, and a new zero-order lattice hamiltonian is obtained, eq. (7.20). To 
understand its physical features, we consider its continuum analog, 

Ho = Hho + Hci (1.16) 

Ho = Ik,;i cok a;tt(k) a;^(k) + (2V)-l Ik P-k^ vk Pk^ (1.17) 

where ^ is a two-valued polarization index, |j, is a hadronic mass, vk = 
k^uk-^, and pk^ and uk are the fourier transforms of the color charge 
density p(x)a and the static color-Coulomb potential defined in (1.3) and 
(1.11a) respectively. Here Hho is a harmonic oscillator hamiltonian with 
frequencies cOk = [k2+|i4uj^]l/2^ ^nd Hci is a static color-Coulomb 
interaction hamiltonian. The color charge density is given by p = 
Pdyn + Pext, where pdyn is the color-charge density of gluons (and of 
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dynamical quarks if any) and pext may be the color-charge density of, say, 
a pair of external quarks. 

The quantity Vk or its fourier transform v(x) has a clear intuitive 
meaning as a color-Coulomb interaction energy between color-charges. 
However v(x) differs from the gauge-invariant energy E(x) of a pair of 
external quarks at a separation x. The latter is the ground-state energy 
Ho^O = E(x)^o, of the system with (approximate) hamiltonian Hq. It also 
appears in the rectangular Wilson loop of dimension T by Ixl, which in the 
limit of large T has the value W = exp[ - E(x)T]. The two differ, 
v(x) E(x), because the external quarks distort the ground-state gluon 
wave-function ^o(At9» by a QCD analog of vacuum polarization, and 
moreover E(x) < v(x), as is shown in Appendix E. 

Because the second term in (1.17) is quartic in the fields, the ground 
state wave function of Hq cannot be calculated exactly. This prevents an 
exact determination of the kernel of the integral equation (6.22) that 
determines the static color-Coulomb potential u(k). In sect. 8, an 
approximate kernel is used in this equation, with solution w(k) = u(0)(k). 
The approximate interaction energy v(0)(k) = k2w2(k) may be written 
v(0)(k) = gc^(k)/k2, where gc(k) is a renormalized running coupling 
constant that characterizes the static color-Coulomb interaction (and the 
approximation made). It is not the gauge-invariant running coupling 
constant g(k) of the renormalization-group, gc(k) g(k), because 
v(x) 4^ E(x), as explained in the preceding paragraph. The asymptotic 
form at high momentum of gc(k) is calculated in sect. 8, and is compared 
in Appendix D with the renormalization group expression for g(k). The 
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two quantities have the same analytic form, but with different leading 
coefficients be ^ brg. Numerically be gives the bulk of the contribution 
to brg, and it is verified that the coefficient be agrees with contribution to 
brg that comes from the conventional perturbative evaluation of <M-1>, as it 
should. At the end of Sect 8, the horizon condition is solved in the same 
approximation for the corresponding unrenormalized coupling constant 
go(k). 

In sect. 9, the asymptotic form at low momentum of the approximate 
static color-Coulomb potential w(k) = u(0)(k) is determined, and thus also 
of the approximate interaction energy v(0)(k) = k^w^(k), which is found 
to have the asymptotic form at low momentum given by v(0)(k) ~ l/kl4/3. 
Such a highly singular infrared behavior implies confinement of color- 
charge in the approximate theory described by Hq and v = v(0). For the 
low-momentum limit of the fourier transformed color-charge density is the 
total color charge p^(k = 0) = Q^, and consequently, by (1.17), there is an 
infrared contribution to the energy of any state given by, 

(^, Ho^) ~ [2(2ji)3]-ljd3k v(k) (W, Q2^) . (1.18) 

This diverges if v(k) is as singular as l/k^, for all states ^ with non-zero 
color charge. The divergence cannot be canceled, because the remaining 
terms in (1.17) are positive. Thus the low momentum behavior found for 
v(0)(k) implies that all states with non-zero total color charge have infinite 
energy. This is a universal mechanism which confines all color charge in 
all representations whether carried by dynamical or external particles. 
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In sect. 10 we present our conclusions. We discuss to what extent the 
confining behavior just described can be expected to survive the 
approximations made, and the relation of the l/kl4/3 behavior of v(k) to a 
linearly rising potential E(x) ~ Klxl. 

2. Restriction of the Configuration Space of a Hamiltonian 
System with a Large Number of Degrees of Freedom. 

The Feynman-Kac formula expresses a quantum mechanical 
expectation value as the expectation value of a fictitious statistical 
mechanical ensemble. If H is the hamiltonian of the quantum mechanical 
system, and Z is the partition function of the statistical mechanical 
ensemble with transfer matrix T, then they are related by 

li exp( - pEi ) = tr[exp(-pH)] = limM^ootrXM, (2.1) 

T = exp(-8H) 1 - £H, for £ = M-lp (2.2) 

Z = trXM. (2.3) 

This correspondence goes both ways: given H we may find T and Z, and 
given T or Z, we may find H. 

In statistical mechanics there are well known equivalences between 
different ensembles such as the microcanonical and canonical. For systems 
with a large number of degrees of freedom, this translates back into an 
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equivalence between different kinds quantum mechanical system. In some 
cases, this allows one to replace a system with a boundary by a much 
simpler system with a modified hamiltonian and no boundary. 

Consider a system with N of degrees of freedom, labeled by an index 
k. In hamiltonian lattice gauge theory, the large-N limit is the infinite- 
volume limit. As a concrete example, the reader may suppose that the 
hamiltonian is quadratic 

H = 2-lIk=lN ( - a2/aqk2 + Ek2 qk2) (2.4) 
where k is a discrete 3 -momentum, and Ek may be, for instance, 

Ek2 = lk|2, (2.5a) 

or 

Ek2 = Ak, (2.5b) 

where -Ak are the eigenvalues of the 3 -dimensional lattice laplacian, or 

Ek = 0. (2.5c) 

However it will be seen that our results hold quite generally for a 
hamiltonian which is bounded below. 
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We shall be interested in the consequences of the restriction of gauge 
field configurations to a fundamental modular region A. As an example, 
the reader may suppose that the fundamental modular region is an ellipsoid 
defined by 

G(q) = Ik=lN qk2/ck2 - N < 0, (2.6) 

where the Ck are all of order unity, but again our results will be seen to 
hold very generally. The wave function ^(q) takes values inside this 
ellipsoid and is required to vanish on the boundary, 

^(q) = for G(q) = 0. (2.7) 

For finite N it is difficult in general to find the eigenvalues of H with this 
ellipsoidal boundary condition, even for E(k) = 0. 

This quantum mechanical problem corresponds to the partition 
function 

Za = Ia nk,t dqk,t exp[ - S(q) ] (2.8) 

S(q) = 2-1 Ik,t [ e-1 (qk,t+l - qk,t)2 + 8 Ek2 qk,!^ ], (2.9) 

where t = 1,... M labels "time" slice. The region of integration is limited 
by the M constraints 

G(qt) = Ik qk,t2/ck2 - N < 0, (2.10) 
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t = 1,... M, as indicated by the subscript A on the integral (2.8). 
To evaluate Za we write 

1a nk,t dqk,t = 1 nk,t dqk,t Dt e[ - G(qt) ], (2.11) 
and we represent the step-function by its fourier expansion 
e[ - G(qt)] = (27i)-l J-oo°° dyt (8+iyt)-l exp[ - iyt G(qt) ]. 

(2.12) 

The integrand is an analytic function of the complex variables zt = xt+iyt. 
For each t, we deform the contour of integration into the right-hand half 
plane, xt > 0, which gives 

Za = j rit [(27ri)-ldzt] exp[ W(z) - It Inzt ], (2. 13) 

where z= (zi,...zm). The path of integration of zt runs parallel to the 
imaginary axis in the right-hand half plane, and W(z) is defined by 

exp[ W(z) ] ^ J nk,t dqk,t exp [- S(q) - It zt G(qt) ] . (2. 14) 

We look for a multi-dimensional saddle point at a critical value, zt = at 
which is real and positive for each t, at > 0. The critical value satisfies 

aw(a)/aat - afi = - <G(qt)>a - afi = 0, (2.15) 
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where cx = (cxi,...cxm), and the expectation value is calculated in the 
ensemble defined by (2.14) with zt = at. The at are all equal as = at. 
The inverse "width" of the multi-dimensional saddle is described by the 
matrix 

32w(a)/3as 3at + 6s,tat-2 



in which the covariance matrix of the statistical mechanical observables 
G(qt) appears. The observable G(qt) is a bulk quantity of the order of the 
number of degrees of freedom N, so <G(qt))a is of order N, and its 
covariance is of the order 5s,tN. Consequently the width of the saddle is of 
order N-l/^ in each variable, and a saddle-point evaluation becomes exact in 
the limit of large N. 

We now prove that (2.15) has exactly one solution with real positive 
a, where a is the common value of the at. As the at are all equal, we may 
replace these M equations by the one equation 



= <G(qs) G(qt)>a - <G(qs)>a <G(qt)>a + 6s,tat-2 



(2.16) 



F(a) = M-l<B>a + a-1 = 



(2.17) 



where 



B = It G(qt) 



(2.18) 
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and M is the number of time slices. We have 



F'(a) = - { M-l [<B2>cc - <B>a2] + a-2 } < 0, 



(2.19) 



so F(cx) is monotonically decreasing. Moreover G(qt) and hence B is 
bounded below, so F(0) = +cx3 by (2.17). At a = +oo, the partition 
function is dominated by the lower bound of H(qt ) which by (2.10) is -N, 
so F(+oo) = - NM < 0. Thus there is precisely one positive solution a. 

We conclude that in the large-N limit, the saddle point 
approximation becomes exact. This gives 

limN^ooZA = limN^ooZa, (2.20) 

where Za is defined by 



where the expectation value is calculated in the ensemble Za. Moreover 
the last equation has exactly one positive real solution. 



Za = I nk,t dqk,t exp [- S(q) - a It G(qt) ], 



(2.21) 



provided that a has the critical value determined by 



F(a) = <G(qt)>a + a-l = 0, 



(2.22a) 
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Remarks : 1 . The critical value a is of order unity, whereas G(qt) is a bulk 
quantity of order N. Consequently in the large-N limit, a may be 
neglected in eq. (2.22a) which simplifies to 



limN^ooN-l<G(qt)>a = 0.. 



(2.22b) 



A further property of a bulk quantity is that its fluctuations are negligible, 
of relative order N-l/^^ so the probability distribution of N-lG(qt) gets 
concentrated at its mean value. 



Consequently when (2.22b) holds, the probability distribution of G(qt) gets 
concentrated at the boundary of the fundamental modular region, at 
G(qt) = in the large-N limit. This is slightly counter-intuitive because 
the probability distribution of the variables qk,t remains peaked at 
qk,t = 0, and G(0) = - N. The paradox is easily resolved however, 
because it is the fluctuation of the qk,t (the mean value of qk,t^) that 
contributes to the mean value of G(qt). We call G the "horizon function", 
and eq. (2.22b) the "horizon condition", where "horizon" refers to the 
boundary of the fundamental modular region. 

2. The result (2.20) is an exact analog of the relation between the micro- 
canonical and the canonical ensemble. 



N-l[G(qt)-<G(qt)>a]-0(N-l/2) 



(2.23a) 



which we write 



G(qt) - <G(qt)>a. 



(2.23b) 
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3. We emphasize that only very general properties of H and G have been 
used in the above saddle-point evaluation. The above results hold for a 
wide variety of systems with boundary conditions, and is not at all 
restricted to the quadratic hamiltonian which is presented as an illustrative 
example only. 

In the proof above that the equation F(a) = possesses a real 
positive solution ttcr > 0, it was implicitly assumed that F((x) is analytic. 
However F(a) may develop a non-analyticity in the large-N limit, and ttcr 
may approach zero, as illustrated in fig. 1. [Indeed, in replacing (2.22a) 
by (2.22b), we dropped the term OC"! which was used to show that 
F(0) > in the proof that the equation F(a) = does in fact possess a 
solution for positive a.] There are two possibilities: 

case (a) <G(qt))o < 0. (2.24a) 

(We write <G(qt)>0 for <G(qt)>a at oc = 0+.) In the large-N limit, the 
fluctuations of G(qt) are negligible, so in the ensemble characterized by 
a = 0+, the probability distribution of G(qt) gets concentrated on the 
surface G(qt) = (G(qt))o ^ that is contained within the fundamental 
modular region. In this case the restriction to the fundamental modular 
region is vacuous, and we have ttcr = 0. 

case (b) <G(qt)>o > 0. (2.24b) 

In this case, the equation (G(qt)>a = holds for some cx = cxcr > 0. The 
fluctuations of G(qt) are negligible in the large-N limit, so in the 
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distribution characterized by a = ttcr > 0, the probability distribution of 
G(qt) gets entirely concentrated on the surface defined by G(qt) = 0, 
namely on the boundary of the fundamental modular region. 

These two cases correspond to different phases in the large-N limit, 
for which <G(qt)) may serve as an order parameter. We call case (a), in 
which the restriction to the fundamental modular region is vacuous, the 
perturbative phase , and we call case (b), in which the probability 
distribution of G(qt) gets concentrated on the boundary of the fundamental 
modular region, the horizon phase . We shall find that in QCD, the horizon 
phase corresponds to the phase of confined gluons. 

In the horizon phase, the entropy is stronger than the energy and 
pushes the system out so that the probability distribution of the horizon 
function is concentrated on the boundary of the fundamental modular 
region. On the other hand, in the perturbative phase, the energy is 
stronger so the probability that G(qt) attains its value on the boundary of 
the fundamental modular region is negligible. 

The equivalence between the quantum mechanical and the statistical 
mechanical system, is obtained in the limit in which the parameter 
8 = M-l(3 approaches zero, where 8 is the parameter that defines the 
length of the "time" interval in (2.9). For small 8, the solution to the 
horizon condition gives a value of a that is of order 8. To see this, observe 
that the canonical ensemble Za, defined by (2.21) with 

a = 8 Y (2.25) 
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corresponds to the quantum mechanical system with effective hamiltonian 



In the large-N limit, the parameter y is determined by the "horizon 
condition" (2.22b) which (in the zero-temperature limit) corresponds to the 
condition 



where is the ground state wave function of the hamiltonian Heff. This is 
independent of 8, so (2.25) is the correct scaling of a. 

The quantum mechanical system defined by Heff with the horizon 
condition (2.27) contains no boundary, but it is nevertheless equivalent in 
the large-N limit to the original system with hamiltonian H and boundary 
condition (2.7). The new system is much more convenient for calculations, 
just as the canonical ensemble is more convenient than the microcanonical 
ensemble. 

To verify that the quantum-mechanical horizon condition (2.27) has 
at most one solution, we observe that F(y) is monotonic. For by first order 
perturbation theory we have 



Heff=H + 7G. 



(2.26) 



F(Y) ^ (^O(Y), G(q) ^o(Y)) = 0, 



(2.27) 



dF(y)/dy = - 2 (^o, G (Heff - Xo)-^ G ^o) < 0, 



(2.28) 
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where A,o = A,o(y) is the vacuum energy of Heff(Y), and the inverse is taken 
on the subspace orthogonal to ^o- 

From (2.26) and (2.27) and the fact that an eigenvalue is stationary 
with respect to infinitesimal variations of its eigenfunction, we have 

F(Y) = 3?io(Y)/3Y- (2.29) 

The horizon condition (2.27) is thus equivalent to the statement that the 
vacuum energy is stationary at the critical value of y. In quantum field 
theory, r(Y) = Xoiy) is known as the effective action or effective potential. 

For the concrete example defined by (2.4) and (2.6), one has 
Heff = 2-1 Ik [ - a2/aqk2 + COkW ] - y N, (2.30) 



where 



0)k2 = Ek2 + 2yck-2. (2.31) 

This is the hamiltonian of a system of decoupled harmonic oscillators with 
eigenvalues 

?i = Ik (nk + 2-l)a)k - yN (2.32) 
where nk = 0, 1 ... . 
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Equation (2.27) reads 



F(Y) = Ik ck-2 (^0, qk^ ^0) - N = Ik (2 cok ck^)-! - N, (2.33) 
so the horizon condition, F(y) = 0, gives 

Ik [ 2 ck ( ck2Ek2 + 2y )1/2 ]-1 - N = 0, (2.34) 
This equation has a solution for positive y if and only if the condition 

Ik(2ck2Ek)-l>N. (2.35) 

is satisfied. If it does, then the system is in the horizon phase with a 
positive value of y satisfying (2.34). If this condition is not satisfied, then 
<G(q))Y is negative at y = 0, and the system is in the perturbative phase and 
the boundary condition may be ignored. In the large-N limit, (2.31) and 
(2.34) provide a complete solution to the original eigenvalue problem with 
an ellipsoidal boundary condition. 

3. Definition and Elementary Properties of the Minimal Lattice 
Coulomb gauge. 

Consider a finite Euclidean periodic lattice which is cubic in the 
space directions. It will be convenient for taking the continuous-time limit 
of a Euclidean lattice to label Euclidean sites by two indices x and t, where 
X is a 3-vector with integer components that labels a spatial site, and t is 
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integer- valued Euclidean time. Spatial, or horizontal, links of the lattice 
will be labeled by (xy, t) and time-like, or vertical, links by (x, t). Lattice 
configurations U are described by horizontal and vertical link variables 
denoted respectively by Uxy,t and hx,t, and local gauge transformations g 
are described by the site variables gx,t. These variables are elements of the 
SU(n) group, and the reversed link corresponds to the inverse group 
element. The gauge-transform of the configuration (U, h) by a local 
gauge transformation g, is written 

(Ug)xy,t = gx,t"'' Uxy,t gy,t. (3.1a) 

(hg)x,t = gx,t"^ hx,t gy,t+l. (3.1b) 

In order to choose a gauge which makes the horizontal link variables 
as close to unity as possible, we require a measure of the deviation of the 
horizontal link variables from unity. For this purpose we take the quantity 

I(U) = It i:(xy) n-1 Re tr(l - Uxy,t) , (3.2) 

where the sum extends over all time-slices t, and horizontal links (xy) of 
the lattice. It is positive, I(U) > 0, and vanishes, I(U) = 0, if and only if 
Uxy,t = 1 on every horizontal link. The restriction of this quantity to the 
gauge orbit through an arbitrary configuration U is given by the (Morse) 
function 

Fu(g) = I(Ug) = It I(xy) n-1 Re tr(l - gx,tt Uxy,t gy,t ) , (3.3) 
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regarded as a function on the local gauge group G, for fixed U. The 
gauge-choice which makes the horizontal link variables as close to unity as 
possible, is achieved by gauge-transforming with the gauge transformation 
that makes this function a minimum. After this transformation, the 
minimum on each gauge orbit occurs at gx,t = 1- The set A of all 
configurations U with this property constitutes the fundamental modular 
region. 



Degenerate absolute minima occur only on its boundary and are identified 
topologically. For a lattice of finite volume, they form a lower 
dimensional set. It is shown below that this gauge choice falls into the class 
of lattice Coulomb gauges, and we call it the minimal lattice Coulomb 
gauge. 

Because the minimizing function (3.3) is independent of the vertical 
link variables hx,t, and is a sum over uncoupled time-slices, the 
minimization independent on each time slice, and the fundamental modular 
region is the topological product 



where At is the fundamental modular region for spatial configurations Ut 
on the time-slice t. 



A^{U: Fu(l)<Fu(g)forallg}. 



(3.4a) 



A = ®t At 



(3.4b) 



At ^ {Ut: Fut(l) < Fut(gt) for all gt}. 



(3.4c) 
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The minimal Landau gauge in D Euclidean dimensions is defined by the 
same condition that holds on each time slice in the minimal Coulomb gauge 
in D spatial dimensions. Consequently the results of [11] for the minimal 
Landau gauge hold here on each time slice. 

We now discuss the character of the minimum on the time slice t, 
temporarily suppressing the Euclidean time variable t. At any local or 
absolute minimum, the first variation of the minimizing function is 
stationary, which gives the so-called gauge condition, and the second 
variation is non-negative. To obtain the explicit form of these quantities, 
we write the local gauge transformation gx = exp(oi)x), where cOx = cOx^t^ 
is an element of the Lie algebra of the SU(n) group. Here the t^ form an 
anti-hermitian basis of the defining representation of this Lie algebra. 
They satisfy [t^, tb] = fabctc^ and are normalized to tr(tatb) = - 2-l6ab. 
We have, to second order in O), 

Fu(g) = Fu(l)-(2n)-ll(xy) 

X [ tr{ [Uxy - Uxyt ] (COy - (Ox + 2-l[0)x, COy]) } 

+ 2-1 tr{ [Uxy + Uxyt ] (0)y - 0)x)2 } ] 

= Fu(l) + (2n)-ll(xy) 

X [Axy'^ (0)yC - COx'^ + 2-lfabco)xa(Oyb) 
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+ 2-1 (0)ya - 0)x^) Gxy^c (coyc . o)xC) ]. (3.5) 



Here we have introduced the real horizontal link variables 

AxyC = - tr[ (Uxy - Uxyt) tc ], (3.6) 

which approach the classical connection in the continuum limit, and 

Gxy^b = . 2-1 tr[ (ta tb + tb ta ) ( Uxy + Uxyt ) ] . (3.7) 

We shall have occasion to also consider lattices with free boundary 
conditions. For this reason it is useful to introduce a notation which is 
inherently geometrical and holds for an arbitrary lattice which may be 
periodic or which has a boundary. For any site variable (px, the link 
variables 

(V(p)xy = (py - (px (3.8) 

(a(p)xy ^ (1/2) (cpx + (PyX (3.9) 

are defined for all links (xy), including those for which x or y may lie on a 
face, edge or vertex. (The lattice operation V, which maps site variables to 
link variables, is obviously a lattice analog of the corresponding continuum 
quantity, called d in Cartan's notation.) As convenient, we also label links 
by (xy) = (x, i), for y = x + ei, where ei points in the positive 
i-direction, so 
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Ai(x) = Ax,y 



(3.10) 



Vi(p(x) =(V(p)xy 



ai(p(x) =(a(p)xy 



(3.11) 



(These expressions are undefined for sites x where y = x + ei is not a site 
of the lattice.) Given two link variables that are group elements Vxy and 
Wxy, we may form the product link variable Pxy = VxyWxy, etc. 

The natural inner product on link variables provides a definition of 
the operator V* which is the dual to V, 



(V*A, (p) = (A, Vcp) = L(xy) Axy (V(p)xy = Ix (V*A)x (px, (3.12a) 



where the sums extend over all links and sites of the lattice. Whereas V 
maps site variables cpx into link variables (V(p)xy, its dual V* maps link 
variables Axy into site variables (V*A)x. We shall also use the notation V- 
for the lattice divergence, V- = -V* 



The lattice divergence is defined for all sites x, and represents the sum of 
all link variables leaving the site x. (It is not formed simply of differences 
of link variables when x is a boundary point of the lattice.) Similarly the 
site variable (a A)x is defined by 



(a- A, (p) = (A, a(p) = L(xy) Axy (a(p)xy = Lx (a-A)x (px. (3.13) 



(V-A)x ^ - (V*A)x . 



(3.12b) 



30 



With the help of these definitions, we rewrite the expansion (3.5), 
which is correct to second order, in the form 



Fu(g) = Fu(l) + (2n)-l [ (A, Vco) + (Vco, D(U)(o) ]. (3.14) 
Here [D(U)co]x,i is the link variable defined by 



It is the infinitesimal change in the link variable Ai^(x) under an 
infinitesimal gauge transformation cOx. By analogy with the continuum 
theory, we call the operator D(U) that maps site variables into link 
variables "the lattice gauge-covariant derivative". We also define the 
lattice Faddeev-Popov matrix M(U) by 



At a minimum, the first variation of the minimizing function is 
stationary, - (A, Vco) = (V-A, co) = 0, for all co, so configurations in 
the fundamental modular region A satisfy the lattice transversality or 
lattice Coulomb-gauge condition 



Diaccoc(x) = Giac(x)Via)C(x) + fabcAib(x)aiCOC(x) . 



(3.15) 



(CO, M(U)(p) = (Vco, D(U)(p). 



(3.16) 



[V-A(U)]x = 0. 



(3.17) 
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At a stationary point of the minimizing function, the matrix of 
second derivatives is symmetric, so when A is transverse, the Faddeev- 
Popov matrix (3.16) is symmetric, 

M(U) = V*D(U) = D*(U)V (for V-A = 0), (3.18a) 

which we write alternatively as 

M(U) = - V-D(U) = - D(U)-V (for V-A = 0). (3.18b) 

Here D*(U) = - D(U)- is dual to D(U), and maps link variables onto site 
variables, and we call D(U)- "the lattice gauge-covariant divergence". The 
Faddeev-Popov matrix is non-negative at a relative or absolute minimum. 
The two conditions together define the Gribov region Q., 

Q. = {U: V A(U) = and M(U) > }. (3.19) 

The set A of absolute minima is contained in the set of relative mimima, 

A cz a. (3.20) 

When properly restricted to A, the Faddeev-Popov formula holds on 
each time-slice 

1 = J dg n 6A[V A(Ug)] detM(Ug). (3.21) 

where 
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5A[V-A(Ug)] = 5[V-A(Ug)] XA(Ug). 



(3.22) 



Here Xa(U) is defined only for configurations U that are transverse, and is 
given by %a(U) = 1 for U in A, and %a(U) = otherwise. By virtue of 
this identity, Wilson's partition function (1.15) may be written in the 
minimal Coulomb gauge 

Z = 1 n dUx,i,t dhx,t 6A[V-A(Ut)] detM(Ut) exp[- Sw(U, h) ]. 

(3.23) 

Here Ut refers to the configuration on the horizontal links of the time- 
slice t, and Sw(U, h) is Wilson's Euclidean lattice action. 

(As discussed in [13], the last expression is symbolic to the extent 
that, for simplicity, we have ignored the fact that on each time slice the 
gauge condition leaves a global or x-independent gauge transformation 
arbitrary. Consequently, on a lattice of spatial volume V, only V-1 gauge 
conditions V A(U)x = are linearly independent, and moreover M(U) has 
a trivial nullspace consisting of constant (x-independent) eigenvectors. The 
determinant detM(U) is defined on the orthogonal space.) 

4. Lattice Coulomb Hamiltonian 

The Euclidean action may be decomposed into parts that lie within a 
given time- slice or link adjacent time- slices. 
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Sw(U, h) = It St+i,t 



(4.1a) 



St+l,t = Pe Ix,i Sx,i(Ut, Ut+l, ht) + Pm Ix,ij Sx,ij(Ut+l) • (4. lb) 



where the variables Ut, and ht live on horizontal and vertical links. The 
sums extend over all vertical and horizontal plaquettes respectively and Sx,i 
and Sx,ij are the corresponding one-plaquette actions. For different 
vertical and horizontal intervals, the coefficients have different values 
given by 



which come from dividing each plaquette action by the square of the area 
of the plaquette, namely (aoa)2 and (a2)2 for vertical and horizontal 
plaquettes respectively, and multiplying by the 4-volume element aoa^. 
(For purposes of renormalization, one may prefer to use different coupling 
constants gt and gs on the vertical and horizontal plaquettes, to avoid 
renormalizing the speed of light [16]. However we shall ignore this point 
here.) 

Only the variables ht appear in St+i,t, but the variables ht+i do not. 
This is the lattice analog of the of absence of 3Ao/3t in the Lagrangian of 
the continuum theory. In the continuum theory, the variable Aq is not a 
true dynamical variable, and variation with respect to Aq gives the 
constraint known as Gauss's law. In this section, the integration over the 
variables hx,t that live on vertical links will be done explicitly in the limit 
of small ao, and will give the lattice form of Gauss's law. 



Pe = 2na/(go2 ao) 



Pm = 2nao/(go^ a). 



(4.2) 
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To find the Hamiltonian, we express the partition function in terms 
of the transfer matrix T by 



Z = tr(TM), 



(4.3) 



where T is related to the Hamiltonian operator Hcoul by 



T = exp( - aoHcoul) = 1 - aoHcoul + OCao^). 



(4.4) 



We now choose the variables Ax,i,t*^ defined in eq. (3.6) as coordinates for 
Ux,i,t on each coordinate patch of the SU(n) group manifold and on each 
time-slice t, so the gauge condition V A = is linear in the coordinates. 
The discrete variable that labels coordinate patches will be suppressed. In 
terms of these variables we have 



where l\|/(Ax,i,t)l is the Haar measure on the SU(n) group in terms of the 
coordinates Ax,i,t [see (A.31)]. As indicated, the integration over dA^^ = 
n^dAtr.cc is restricted to the fundamental modular region A, where, for 
each time-slice t, the Atr,a are a set of independent coordinates that 
parametrize the transverse components of Ax,i,t. They will be chosen to be 



1 11x4 dUx,i,tnx6A[(V-A)x,t] = 



J 11x4 dAx4,t l¥(Ax4,t)l nx6A[(V-A)x,t] 



= Ja dAtr 0x4 l¥(Ax4,t^OI , 



(4.5) 
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the fourier coefficients of Ax,i,t. (The Atrx^^t are not independent variables 
on each Hnk.) 



From the partition function (3.23) and (4.1), one reads off the 
transfer matrix T, 

(TO)(Atr) = N Ja dBtr detM(AttO Oxdgx 

X exp[- St+i,t ] ^(Btr), (4.6) 
where N is a normalization constant, 

St+l,t ^ Pe Ix,i Sx,i(Atr, Btr, g) + Ix,ij Sx,ij(Atr). (4.7) 
and we have written 

Ax,it^ = Ax,i,t+lti-; Bx,iti- = Ax,i,ttr; ^ hx,t • (4.8) 

The horizontal one-plaquette action Sx,ij = Sx,ij(Atr) depends only 
on the horizontal variables within a single time slice. By comparison with 
(4.4), one sees that the horizontal plaquettes contribute a potential energy 
term to the hamiltonian, which, with j3m = 2nao(go^a)-l, is given by 

V(Atr) ^ 2n(go2a)-l Ix,ij Sx,ij(AtO. (4.9) 

We write 
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Hcoul = K + V(Atr 

or 

T = exp[ - aoV(Ati-) ] exp[ - ao K ] . 
The remaining part of the transfer matrix 

T' = exp(- aoK) = 1 - aoK + O(ao2) 
acts according to 

(T'VF)(Atr) = N Ia dBtr nx,i \\\f(B^,i^^)\ detM(Btr) nxdg 
X exp[- 8 -1 Ix,i Sx,i ] ^(Btr) . 
Here the vertical one-plaquette action Sx,i is given by 

Sx,i = 1 - n-1 Re tr [Ut(Ax,itr) gxt U(Bx,itr) gy ], 
where y = x + ei, and 

8 = go^ ao (2na)-l 
is a small quantity that vanishes with ao. 

37 



The evaluation of K, which turns out to be the kinetic energy 
operator, is done in Appendix A. The most transparent way to express K 
and the Coulomb Hamiltonian Hcoul is as quadratic forms. Let O(Atr) and 
^(Atr) be wave functions, where we have restored the superscript tr for 
transverse. The inner product is given by 

(<D, ^) = Ja dAtr (y(Atr) <D*(Atr) ^(Atr), (4. 1 5) 

where a(Atr) = det\|/(Atr) detM(Atr), and M(Atr) = D*(Atr)V. The sub- 
script A indicates that the integration is restricted to the fundamental 
modular region. With kinetic energy operator (A. 54) and (A. 55), and 
potential energy (4.9), the quadratic form corresponding to the Coulomb 
hamiltonian is given by 

(O, Hcoul^) = Ja dAtr (y(Atr) { [ go2(2a)-l] Ix,i (£x,i^^)*(£x,i^^) 

+ 2n(go2 a)-10* [ Lx,ij Sx,ij ] ^ }, (4.16) 

where <I> = <I>(Atr), ^ = ^(Atr), and Sx,ij = Sx,ij(Atr). Here £x,i^ is the 
lattice color-electric field operator defined by 

^x,i,a ^ JaKAti-x,i) Ex,i,b (4.17) 

where Ja'^(Ax,i) is given in (A.56) or (A.57), and Ex,i,b acts according to 

= i [ a^/aAtr . V (D*V)-1 D*a^/aAtr ]. (4.18) 
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The electric field operators satisfy the lattice non-Abelian Gauss's law 
constraint, 

Z)*£ = D*E = 0. (4.19) 
The electric field operator may be written 

E = Etr - V(p, (4.20) 

where 

Etr = i a/a Atr, (4.21) 
and cp is the non-Abelian lattice Coulomb potential operator. 

9 = i (D*V)-1 D*a/aAtr. (4.22) 
This completes the calculation of the lattice Coulomb hamiltonian Hcoul- 

Remarks : 1. In case color charge is carried by other fields such as the 
quark field, then expression (4.22) for the color-Coulomb operator is 
modified to 

(p = (D*V)-l ( i D*a/aAtr + pqu), (4.23) 

where Pqu is the quark color charge density operator. Gauss's law is 
modified to 
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- D*E = D E = pqu. 



(4.24) 



Gluons contribute to the color-Coulomb operator by the charge density 
operator iD*a/aAtr = -iD-d/dA^^ 

2. We give a more explicit expression for the case of SU(2). In this case 
we have on a given link 

U(A) = r| [ 1 - (A/2)2 ] 1/2 + taAa, (4.25) 

where we have written for At^x^i^, and r| = ±1 is a sign that labels 
coordinate patches on the SU(2) group manifold. From UtdU = t^^\dA^, 
we obtain for the Maurer-Cartan form for the SU(2) group in these 
coordinates, 

\|/(A)ac = r| [ 1 - (A/2)2 ]l/2 6a,c 

+ 4-1 r| [ 1 - (A/2)2 ]-l/2 AaAc - 2-lfa,b,c Ab, (4.26) 

where fa,b,c ^re the structure constants of the SU(2) group. From 
\|/acJab = 6cb, we obtain 

J^c = n [ 1 - (A/2)2 ]l/2 - 2-lfa,b,c Ab, (4.27) 

(J*J)b,c = j^b j^c = §b,c . 4-lAb AC, (4.28) 

and the Haar measure on a single link is 



40 



detv|/ = 11- (A/2)2 1-1/2 



(4.29) 



Equation (4.28) for J*J on a given link yields for the sum in (4.16), 
Ix,i (£x,i^^)*(^x,i^^) = Ix,i [ 6b,c . 4-lAtrx,ib Atr^^c ] 

X (Ex,ib<J))*(Ex,i'^^), (4.30) 

where Ex,ib is defined in (4.18). 

3. We have chosen a lattice notation which agrees with the continuum 
notation for corresponding quantities. To obtain the continuum expression 
for the kinetic energy operator, it is sufficient to set J=l, and 
correspondingly \|/ = J*-l. Our expression then agrees with the continuum 
hamiltonian of Christ and Lee [2]. In the continuum, the gauge-co variant 
derivative is given by Di^'C = 6^''^Vi + f^bcAtr^b^ so the contribution of the 

gluons to the charge density operator which appears in the Coulomb 
potential is i (D'^d/dA^^)^ = - i (D-d/dA^T^)^ = - i fabc At^ib d/dA^^f. 

4. The kinetic energy operator K is not the Laplace-Beltrami operator of 
the metric G, which is the operator that would be obtained if o = 
det\|/ detM were replaced by det^/^G = a dQt-^^^(D*D), in our expression 
for K. Of course the measure of integration in (4.15) and (4.16) could be 
trivially changed from o into detQl/^^ by rescaling the wave functions by 
the factor det^^^(D*D). But then the operator K would be changed into the 
Laplace-Beltrami operator of the metric G plus a non-local potential term 
involving (Z)*Z))-1. The interesting point is that it is possible to write Hcoul 
without any non-locality associated with det(Z)*Z)) or (/)*/))-!, even though 
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these quantities appear in intermediate stages of the calculation. The only 
non-locality Hcoul is associated with the inverse Faddeev-Popov operator 
M-l(Atr) = D*(Atr)V. 

5. In the classical Yang-Mills action, 

S = 2-lJd4x F^iva(aj^Ava - 3vA^a + ggfabcA^bAvC) - 

variation with respect to Aq produces Gauss's law, D-E = 0. This is the 
classical analog of our lattice procedure of integrating out the variables that 
live on the vertical links, which produced the lattice version of Gauss's law. 

6. We would have obtained the same lattice Coulomb hamiltonian if we 
had followed an alternative procedure which is the lattice analog of the 
continuum method of Christ and Lee [2]. These authors first calculated the 
hamiltonian Hw in the "Aq = 0" or Weyl gauge, where the continuum 
coordinates Ai(x) are Cartesian, and where there is a residual invariance 
under time-independent gauge transformations. They then obtained the 
continuum Coulomb hamiltonian Hcoul by making the non-linear change of 
variables Ai = g-lAi^oulg + g-lVig, and then letting Hw act on gauge- 
invariant wave-functions. (The use of Cartesian coordinates as an 
intermediate step was advantageous for these authors whose purpose was to 
quantize a classical theory. In the present article, Cartesian coordinates are 
not needed, because we are gauge fixing the theory which is already 
quantized and regularized by Wilson.) The lattice analog of the "Aq = 0" 
gauge is the "Uo = 1" gauge. The lattice hamiltonian in the "Uq = 1" 
gauge is the Kogut-Susskind hamiltonian [15], which is of the form (4.16) 
and (4.17), with the following modifications: the variables Ax,i^ 
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independent on each link; a (A) = det\|/(A); the color-electric field 
operator Ex,i,b is given simply by 



Ex,i,b = i 3/3Ax,i (4.31) 

instead of (4.18), and there is no restriction to a fundamental modular 
region. However in the "Uq = 1" gauge, there is a residual gauge- 
invariance under time-independent gauge transformations, and wave- 
functions are required to satisfy Gauss's law as a subsidiary condition to 
satisfy gauge invariance. (The left-hand side of Gauss's law is the 
generator of time-independent gauge transformations.) This is an 
independent assumption in the "Uq = 1" gauge, whereas in the Coulomb 
gauge, we have derived gauge-invariance in the form of Gauss's law, by 
integrating over the vertical links. The last sentence contains a slight over- 
simplification. The quantum mechanical partition function Z = 
tr exp(-PH) is a trace, and this in fact corresponds to periodic boundary 
conditions in Euclidean time. Strictly speaking, the "Uq = 1 " gauge does 
not exist with periodic boundary conditions, because it assigns to unity the 
value of Wilson loops that close by periodicity. A correct procedure is to 
make all vertical links unity except those at a fixed Euclidean time, and one 
may show that integration over these last vertical links indeed effects the 
projection onto gauge-invariant states. However we have preferred not to 
first pass to the Uq = 1 gauge in order to maintain the connection with 
"optimal gauge fixing", and because our statistical mechanical arguments 
do not allow a special role for a fixed Euclidean time slice. 
7. For the fundamental modular region A to be truly the quotient- space of 
configurations U modulo gauge transformations g, it is necessary to 
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identify some points on the boundary 3A which are gauge-transforms of 
each other. This means that the wave-functions are required to assume the 
same value at these points. However, because the Coulomb gauge may be 
regarded as the "Uq = 1" gauge, with the subsidiary condition that the 
wave-functions be gauge invariant, the wave functions that are 
eigenfunctions of the Coulomb hamiltonian satisfy this requirement 
automatically [5]. Thus we need not be concerned to identify these points 
explicitly, and it is sufficient to cut off the integration over A at the 
boundary 3 A. 

8. The Gauss's law constraint satisfied by Ex,i^, the definition of A^, eq. 
(1.7), and the Heisenberg equations of motion of A = A^r and Rtr 
generated by the hamiltonian Hcoul provide a quantized and regularized 
form of all the classical Yang-Mills equations on a lattice of finite spatial 
extent. 

5. Fundamental Modular Region and Effective Coulomb-gauge 
Hamiltonian 

The Coulomb hamiltonian, (4.16), does not by itself give a complete 
expression of the dynamics of a gauge theory, because the space of wave 
functions must also be restricted to the fundamental modular region A. 
(Alternatively, a signed measure may be used, as discussed in [4] and [5]). 
An explicit expression for its boundary dA is not known for a lattice of 
finite volume. Fortunately, a limiting expression for 3A, valid at large 
volume, was derived for the minimal lattice Landau gauge on a periodic 
lattice in D Euclidean dimensions in sects. 2-4 of [11], which holds also for 
the minimal Coulomb gauge on a given time slice in D spatial dimensions. 
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(The horizon function which was designated H in [11] is designated G in 
the present article.) 



Let U = U(Atr) be ^ transverse configuration on a single time slice 
of a periodic lattice. Let G(Atr) \jq defined by 



where, a is the lattice spacing, and M = D*V = V*D is the Faddeev- 
Popov matrix. Its inverse M-1 is taken on the subspace orthogonal to the 
trivial null space consisting of constant (x-independent) vectors. The lattice 
difference operator V and the lattice gauge covariant derivative 
D = D(At9, and their duals V* and D*(At9 are defined in sect. 3. (The 
sum in the second term is the minimizing function (3.2), apart from an 
additive constant.) Let A be the fundamental modular region, namely the 
set of configurations U(At9 which are absolute minima of the minimizing 
function I(U), eq. (3.2), on each gauge orbit, I(U) < I(Ug). Let A' be the 
subset of the Gribov region Q defined by G(Atr) < q, 



aG ^ Ix,y Dx,iac Dy,iad m-1 



x,c;y,d 



(n2-l) I(xy) n-1 RetrU: 



xy , 



(5.1) 



A' = {Atr; G(Atr) < 0} n Q . 



(5.2) 



Let Za and Za' be the corresponding partition functions 



Za = Ia n dAtr a exp[ - S(U) ] 



(5.3) 
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and 



Za- ^ In n dAtr a e[ - G(U) ] exp[ - S(U) ], 



(5.4) 



where o = det\|/detM, det\|/(Atr) is Haar measure, and S is a generic action. 
Then these partition functions (and the corresponding expectation values) 
are equal in the infinite- volume limit, 



Loosely we may say that A = A' in the infinite-volume limit. 
Remarks : 

1. In the interior of the Gribov region the horizon function G(At^") is 
finite, but approaches infinity as A^^ approaches a generic point on the 
boundary dO.. Proof : By definition of Q., the Faddeev-Popov matrix 
M(Atr) is strictly positive (apart from the trivial null space) for all A^^" in 
Q, and it develops a (non-trivial) null eigenvalue on the boundary of Q.. 

2. For all U(Atr) i^ ^2, the horizon function G(U) is bounded from below 
by its value at the vacuum configuration Uxy = 1, namely 



where V is the spatial volume and D is the number of spatial dimensions. 
(The quantity F = (n2-l)VD is the total number of degrees of freedom of 
link variables on the lattice.) 



limv^ooZA = limv^ooZA'. 



(5.5) 



G(U) > G(l) = - (n2-l)VD, 



(5.6) 
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Proof : The first term in (5.1) is positive for all configurations U in the 
Gribov region, because the Faddeev-Popov operator is positive there. This 
term attains its minimum value, zero, for the vacuum configuration 
Uxy = 1, because Dx,i(U) is the lattice difference operator at U= 1, 
Dx,i(l) = Vi, where Vi is the lattice difference operator, and on a periodic 
lattice 

Lx (Vif)x = 0, (5.7) 

where fx is any site variable. The second term of G(U) is obviously also a 

minimum at Uxy = 1, where it has the stated value. 

3. The horizon function G(Atr) vanishes identically on a finite lattice with 
free boundary conditions for any transverse configuration Atr. Proof : Let 
i and a be a pair of indices with the ranges i = 1,...D, and a = l,...(n2-l), 
and let fi,a be the color-vector valued site variable defined by 

fi,a,xb = Xi 5ab. (5.8) 

This function is well defined on an infinite lattice or a finite lattice with 
free boundary conditions, but it is not well-defined on a finite periodic 
lattice. One has 

[M(Atr)f]i,a = - D(Atr).[Vfi,a] - - D(Atr).(ei Ca), (5.9) 

where Ca is the color-vector with components = 6^^, and ei is the link 
variable with the value 5ij on links (x,j). With this expression for 
[M(AtOf]i,a it is easy to verify that 
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G(Atr) ^ [ (Mfi^a, M-1 Mfi,a) - (fi,a, Mfi,a) ] = 0. (5.10) 

This equation is false on a finite periodic lattice, as shown by (5.6), but it 
holds rigorously and identically on a finite lattice with free boundary 
conditions. To the extent that the boundary of the lattice becomes 
irrelevant in the infinite- volume limit, the equation G(Atr) = q j^ay be 
considered to hold in a formal sense on a periodic lattice of infinite volume 
when the horizon condition V-1<G> = is imposed as a limiting condition 
at large volume. 

4. It might be thought that the introduction of the term yoG into the 
hamiltonian would destroy renormalizability. Renormalizability is not well 
explored in the Coulomb gauge. However it should be noted that the 
analogous procedure in the Landau gauge does not destroy 
renormalizability for any value of yo. In that case, one adds yoG to the 
Euclidean action, and this non-local term is made local by an integral 
representation in terms of auxiliary ghost fields. It has been shown [14] 
for any value of yo, that in a perturbative expansion in the Landau gauge, 
the renormalization constants may be chosen independent of yo, and thus 
have the same value as in theory with yo = 0. The basic reason is the 
preceding remark, which shows that G is, in essence, a boundary term. 

The horizon function G(A^^) is a bulk quantity of the order of the 
lattice volume, so the results of sect. 2 are valid in the infinite-volume 
limit. Consequently the statistical mechanical system described by the 
partition function Za' is equivalent at large spatial volume to the quantum 
mechanical system described by the effective hamiltonian 
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Heff = Hcoul + yoG, 



(5.11) 



where Hcoul is given in (4.16) and G in (5.1). In the formal sense of 
remark 3, we have G ~ at infinite- volume, so in this formal sense 
nothing has been added to HcouL and we have Heff ~ Hcoul- 

To calculate with this hamiltonian one makes the usual rescaling 



Having taken into account the boundary of the fundamental modular 
region, we shall ignore other possible non-perturbative effects, such as the 
existence of different coordinate patches to cover the SU(n) group 
manifold, and the integration on A^'^^ will be extended to infinity. 

It is convenient to rescale the wave function according to 



where a = det(\j/M), which gives the inner product the simple form 



Atr = go AO,tr, 



(5.12) 



(5.13) 



(O, ^) = In dAtrO'* 



(5.14) 



and the hamiltonian becomes 



Heff' = al/2 Heff a-1/2. 



(5.15) 
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The advantage of this transformation is that o = det(V|/M) appears only in 
Heff' in the combination 



O-l ao/aAtr = a[tr ln(\|/M)]/aAtr 



= tr(M-iaM/aAtr) + tr(\|/-ia\|//aAtr), 



(5.16) 



so it is only necessary to calculate the inverse M-1 instead of the 
determinant, detM. The explicit form of the hamiltonian after this 
rescaling is given in Appendix B. 

6. Color-Coulomb Potential 

In the remaining sections A is rescaled according to A = goA(O), and 
d/dA = go"ia/aA(0), and the superscript on A(0) will be dropped. We also 
set the horizontal lattice spacing a to unity. 

We define the 3 -dimensional Faddeev-Popov propagator Cx and its 
fourier transform C(0) by 



Cx-y 6a,b = <M-1 




(6.1) 



Cx-y = V-lIe exp[ie-(x-y)] c(e), 



(6.2) 



50 



where M is the Faddeev-Popov matrix, and we use the same symbol for the 
fourier transform. In Appendix C we show that C(0) is of the form 

C-l(e) = A(0) - qiO) EijO) qjO), (6.3) 



where 



qi(e) =2sin(2-lei) (6.4) 

AO) ^ li qe,i2 . (6.5) 

The second term in (6.3) is a self-energy, so Z is a self-energy with 
momentum factors removed. We also show in Appendix C that horizon 
condition is equivalent to 

[ E(0)ij - 6ij ]TT = 0. (6.6) 

We suppose that the integral for Sij(O) is (infra-red) convergent, as 
will be verified a posteriori, (8.22) below. Cubic or rotational invariance 

gives 

Eij(O) = const. X 6ij, (6.7) 
and the horizon condition (6.6) is equivalent to the statement 

Sij(O) = 6ij. (6.8) 
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From (6.3) this gives the important result, 

C-l(e) = qi(e) [ ZyCO) - ZyCe) ] qj(e). (6.9) 

Thus the restriction to the fundamental modular region expressed by the 
horizon condition implies that the Faddeev-Popov propagator, which 
transmits the static color-Coulomb potential, is highly singular at k = 0. 

We now come to a crucial point. The last equation contradicts the 
usual perturbative expansion, according to which E is of leading order go^, 
whereas the left hand side is of leading order go^. To find a different 
expansion in go that is compatible with the last equation, we use the 
expression for the self-energy derived in Appendix C 

i:e,ij =go2 [(n2-i)v3]-i ie-,e" cos(2-ie'i) cos(2-ie"j) 

X fabc faed < Ae.e',ib Ae".eje (M-l)e',c;e",d > + .... (6.10) 

where the dots represent higher order terms, and the fourier transform of 
M-1 is defined by 

(M-l)e',c;e",d = Ix,y (M-l)x,c;y,d exp[i(- G'-x + 0"-y)], (6.11) 

For 0' = 0", this gives 

V-l(M-l)e',c;e',d = Iz exp( i0'-z) V-1 Ix(M-l)x,c;x+z,d- (6-12) 
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For each z, the quantity V-1 Lx(M-l)x,c;x+z,d, being an average over all 
sites of the lattice, is of the form B/V, where B is a bulk quantity. Thus, as 
discussed at the end of sect. 2, if we separate out its vacuum-expectation or 
mean-field value, the remainder is of order 

V-1 Ix(M-l)x,c;x+z,d = V-1 < Ix(M-l)x,c;x+z,d> + 0(V-l/2) 

V-1 i:x(M-l)x,c;x+z,d = Cz 5c,d + 0(V-l/2) (6.13) 

V-l(M-l)e,c;e,d = Ce 6c,d + 0(V-l/2). (6.14) 

We conclude that in the infinite-volume limit the part of M-1 that is 
diagonal in momentum space may be equated to its mean- value. 

Consider the contribution of this diagonal part of M-1 to the self- 
energy (6.10) 

^e,i,}=so^ [(n2-i)v2]-i ie-,e" cos(2-ie'i) cos(2-ie"j) 

X fabc faec < Ae.e',ib Ae".eje ) Ce' + ... . (6.15) 

This is a disconnected part of I^e,i,j in the sense that it is a product of 
vacuum expectation values. By translational invariance only the term with 
0' = 0" survives 

i:e,ij = n go2 v-1 Le' De-e-y cos(2-i0'i) cos(2-i0'j) Ce- + ... . 

(6.16) 
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where 



De,i,j 6b,e ^ £y exp(-ie-y) < Ax,ib Ax+yj^ > (6.17) 

is the equal-time gluon correlation function in momentum space. We 
substitute this expression into (6.9), and obtain 

Ce-1 = n go2 v-1 le- qi(e) [ De-y - De.e-ij ] qj(e) 

X cos(2-ie'i) cos(2-ie'j) Ce' + ... . (6.18) 

This is an integral equation for Ce. Because of the coefficient go^, it 
blatantly contradicts the usual perturbative expansion, according to which 
C and D are of leading order go^. To solve this equation for Ce, we 
require an expansion in powers of go which is consistent with it. We 
therefore suppose that the Faddeev-Popov propagator Ce is of leading 
order go'^, 

Ce = go-iue + 0(i). (6.19) 

Here ue is independent of go, and we call it "the static color-Coulomb 
potential". By (6.14), M"! acquires a mean-field contribution that is also of 
order go'^, 

v-i(M-i)e,c;e',d = go"^ ue 6e,e' 5c,d + v-iRe,c;e',d, 
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(M-l)x,c;y,d = gO"^ Ux-y 6c,d + Rx,c;y,d. 



(6.20) 



We make the Ansatz that R and all other correlation functions are power 
series in go with non-negative powers. 

With this Ansatz, the terms in (6.18) represented by the dots are of 
higher order in go. Upon equating terms in (6.18) that are of leading 
order in go, we obtain an integral equation for u, 

u-l(e) = n (2ji)-3 Jd30' qi DOij(e') qj { cos(2-le'i) cos(2-l0'j) u(e') 

- cos[2-i(e'+ e)i] cos[2-i(e'+ e)i] u(e'+ e) }. (6.21) 

where qi = qi(6), and we have taken the infinite- volume limit. Here 
D^i,j(6') is the equal-time gluon correlation function calculated to zero- 
order in go. This non-linear integral equation determines the static color- 
Coulomb potential u(0), once D^j j(6) is given. It holds only in a non- 
abelian theory because the coefficient n which appears in front of the 

integral comes from the contraction of two structure constants, 
fabcfaec = nS^e. 

Because of the difference that appears in the integrand, this equation 
remains finite in the critical limit, where it simplifies to 

u-l(q) = n (2jl)-3 Jd3k DO(k) [ q2 - (q.k)2/k2 ] [ u(k) - u(k+q) ] . 

(6.22) 
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The horizon condition fixes go = go(^)- The form (6.8), or 
equivalently, 



i:i,i(0) = 3, 



(6.23) 



is most appropriate because, by (6.3), it assures that Ce is consistent with 
our expansion in every order. This gives, in leading order. 



gO-1 = n [3(271)3]-! 1 d3e li DOi^e) cos2(2-lei) u(e), (6.24) 
where we have used (6.18) and (6.19) for l^i,i(0). 



7. Re-expanded Zero-order Hamiltonian 

We require the kernel DO of the last equation which is given in 
(6.17). To evaluate this expectation-value to zero order in go we need the 
ground- state wave-function ^o of the zero-order hamiltonian Ho. 
However to obtain the zero-order hamiltonian, we must use the new 
expansion, (M-l)x,y = go"^u(x-y) + ... wherever M"! appears in the 
Hamiltonian Heff' = Hcoul' + YoG. It appears in Heff' (B.l) or (B.3) in 
three places. 

First, it appears in the Coulomb contribution to the color-electric 
field operator (4.18), 



E = Etr - Vcp, = i a/aAtr . gQ VM-lp, 



(7.1) 
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where the color-electric charge density operator is defined by 



p = i B*a/aAtr. (7.2) 

Here we have replaced D* by goB* = D* - V*, because V*3/3Ati' = 0. 
To leading order in go, the color-electric field and color-charge density 
operators are given by 

Ex,ia = i a/aAx,ia,tr . Vx,i ly Ux-y py^ + ... (7.3) 

p^a = . fabc 2-1 £i [Aitr,b(x) Eitr,c(x) + Aitr,b(x-ei) Eiti-,c(x-ei) ] + ... 

(7.4) 

where the dots represent terms that are higher order in go. In the new 
expansion, the Coulombic color-electric field operator has a zero-order 
piece, just like Etr. 

Second, M-1 appears in the horizon function yoG, eq. (5.1), which 

we write 

YoG = YO [ gO^ Ix,y Bx,iab By,iac M-lx,b;y,c 

- (n2-l) Ix,i n-1 Re tr Ux,i ] . (7.5) 

Here we have replaced D by goB, because LxV vanishes with periodic 
boundary conditions. We have 
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yoG = yo[-3(n2-l) V 



+ go2 Ix,y f^db faec Axjd Ay4e ai(x) aj(y) M-lx,b;y,c + ... ] (7.6) 
which gives 

yoG = yo [ - 3 (n2-l) V + n go Ix,y Ax,ib Ay,ib aj(y) Ux-y + ... ], 

(7.7) 

where ai(^) is defined in (3.9). The order in go of these terms depends on 
the order in go assigned to yo. We pose 

yo = (2ngo)-l (7.8) 

where |i is of zero-order in go, so that yoG gives a non-trivial contribution 
to Ho. A perturbative expansion will be made in powers of go, at fixed |i. 
Subsequently the horizon condition will be used to fix go(|-i). This is 
similar to the substituting the solution of the renormalization-group 
equation for go into the usual perturbative expansion. 

By this assignment, yoG has the expansion, 

yoG = -3(n2-l) (2 n go )-l |i4 V 

+ 2-l|i4 £^ y ib Ay,ib ai(x) ai(y) Ux-y + ... . (7.9) 
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The constant term of order go'^ is irrelevant for the hamiltonian and will 
be dropped from Hq, although it contributes to the horizon condition 
<G> = 0. 



Finally M"! appears in the term 



tr(M-iaM/aAx,i^) = Iy,z (M-l)y,b;z,c ^^c a(goAd*V)z,y/aAx,i^ 



+ ... 



(7.10) 



in (5.16) that is present in Heff', and comes from absorbing the integration 
measure detM in the wave function. However because of the trace on the 
color indices, (7.10) has no contribution of order go^. Thus to zeroth 
order we may replace Heff' by Heff. 

To write the hamiltonian explicitly, independent variables must be 
chosen so the transversality condition is satisfied identically. This is 
accomplished using the fourier components defined by 



A^\i = (2V-l)l/2 I'e,^ 8e,i^{ Aec^ cos[e-(x + 2-lei)] 



+ Aq^'^ sin[e-(x + 2-lei)], 



(7.11) 



= (2V-l)l/2 I'e,;^ £e/{ Ee^A cos[e-(x + 2-lei)] 



+ Ees,^ sin[e-(x + 2-lei)], 



(7.12) 
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where the color index is temporarily suppressed, 9i are the lattice 
momenta, Gi = 27ini/N, for ni = 1, 2,...N, and the prime on the sum over 
0i indicates that the sum extends over half the Brillouin zone, and the 
polarization vectors 894^ will be specified shortly. To avoid a profusion 
of the superscript tr, the symbols A and E here and below refer only to the 
transverse components, which are the independent dynamical variables. 
These fourier components satisfy the canonical commutation relations 

[EetA, Ae-t''^'] = - [Pe^'^, Ae-f'^'] = i 6e,e' 5t>t' 5^'^' , (7.13) 

where t = c or s. 

We introduce the exponential fourier transforms 

Ab,trx4 = V-1/2 Ze,X ee,i^ Aeb>^ exp[ie-(x + 2-lei)] 

Eb,tr^ i = V-1/2 ^e^^ ee Eeb,X exp[- ie-(x + 2-lei)], (7.14) 

where the color index is restored, the sums now extend over the whole 
Brillouin zone, and 

Ae?i = 2-1/2 (Aec,^ - i Ae^'^); A.e^ = (Ae^)* (7.15) 
Ee?i = 2-1/2 (Eec,;i + iEesA); e.g^ = (Ee^)*. (7.16) 

These quantities satisfy the commutation relations 
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] = i 6e,e' 5b,b' 5x,x'. 



(7.17) 



The transversality conditions on Atr, namely Li [ Ati'i(x) - Ati'i(x-ei) ] 
= 0, where Atri(x) = Atr^ i, and similarly for Etr, are easily satisfied in 
momentum space because 

exp[ie-(x + 2-lej)] - exp[ie-(x - 2-lej)] = iq(ej) exp(ie-x) (7.18) 

where q(6i) = 2sin(2-l0i). It is sufficient to choose the polarization 
vectors to be a pair of orthonormal polarization vectors that are 

orthogonal to q(0i), 



To obtain the hamiltonian to zeroth order in go we write expressions 
(4.16), (5.11), (7.3), (7.4) and (7.9) in terms of the canonical variables just 
defined, with the result 



Hho = 2-110 {E.eb,X Eeb,X 

+ [ Ae + |i4 he^,^ u(e) ] A-g^.k aq^^^ ] , (7.21) 



Ii=l3 q(ei) 894^ = ;i=i,2. 



(7.19) 



Ho = Hho + Hci, 



(7.20) 



where 
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is a harmonic oscillator hamiltonian, and Hci contains the color-Coulombic 
interaction. The first term in (7.21) is the usual zero-order kinetic energy 
of the transverse color-electric field. The second term is a potential energy 
which is the sum of the usual zero-order potential energy of the color- 
magnetic field with coefficient Ae = Li qe^^ plus the zero-order potential 
energy of the horizon function yoG that contains the color-Coulomb 
potential u(9). The matrix Le^.A, is defined by 

Le^'^ ^ li 8e,i^ cos2(2-iei) ee/ ] . (7.22) 

In the critical limit, Lq^^^ = 5^'^. 

The color-Coulomb interaction hamiltonian Hci is the contribution of 
the zero-order Coulombic electric field (7.3), 

Hci ^ 2-1 Ix,i (Vi(px)2 = 2-1 Ix,i (Vi ly Ux-y Pyb)2 

Hci = (2V)-iIe p-eb ve pe^ . (7.23) 

where the color-electric charge px^* is given in (7.4). It represents the 
color-Coulomb interaction energy with potential 

ve = Ae u2e , (7.24) 

The fourier components pe'^ of the color-charge operator (7.4) are given 
by 
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p^a = V-1 Ix exp(ie x) pe^ 



(7.25) 



pea = . fabc Be+e',e'^'^' Ae+e-'^'^ Ee'C>^' , (7.26) 

where 

Be,e'^'^' ^ li ee,i^ cos[2-i(ei-ei')] ee',i^'. (7.27) 

In the critical limit this simplifies to 69^9'^'^' = 89^ £9'^ . (In 
accordance with our convention, A9 and E9 are components of the 
transverse fields only.) The zero-momentum component of is the 
total color-charge operator 

Qa = p9=o^ = - f^bc ^0 X^hX E9C,?i (7.28) 

that satisfies the commutation relations of the Lie algebra of SU(n), 

[Qa Qb] = i fabc Qc. (7.29) 

Upon diagonalizing the 2 by 2 matrix (7.22) one obtains harmonic 
oscillator frequencies 0)9,x,, and eigenpolarizations 894^. We shall not do 
this explicitly, but note only that in the critical limit the harmonic oscillator 
frequencies are given by 

0)(k)= [k2 + |i4u(k)]l/2. (7.30) 
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If the gluons were free particles, this dispersion relation would contradict 
special relativity (unless u(k) is a constant), but if the gluons are confined 
by the color-Coulomb force, then there is no obvious contradiction. 

With the introduction of the creation and annihilation operators, 

ae^.^ = (2coe,>.)-i/2 ( coe,>. Aq^-^ - i E.eb,^ ) 

(aeb>^)t = (2coe,;i)i/2 ( coe,;i A-q^-^ + i Eeb>^ ), (7.31) 

which obey the commutation relations 

[ (aeb'^)t, ae-b''^'] = 6e,e' 5b>b' (7.32) 

the harmonic oscillator hamiltonian takes the standard form 

Hho = le,^ coe,^ [(aeb,^)t aeb,^ + 2-I]. (7.33) 

If quarks are present, the zero-order hamiltonian becomes 

Ho,tot = Ho + Ho,qu, (7.34) 

where Ho,qu is the free quark hamiltonian. Here Hq is as above, but with 
color-charge density defined by 

p = Pgl + pqu, (7.35) 
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where pgi is given in (7.26) and pqu is the color-charge density of the 
quarks. The color-Coulomb interaction potential v(9) couples at low 
momentum to the total color charge of all colored particles. If the 
interaction potential ve is sufficiently strong at zero momentum then the 
energy of all states with non-zero color charge is infinite. 

8. Determination of the Short Range Static Color-Coulomb 
Field 

The static color-Coulomb potential, u(k), is the solution of the 
integral equation (6.22) whose kernel is the equal-time gluon correlation 
function D^, eq. (6.17), evaluated in the ground state of 
Ho = Hho + Hci. However cannot be calculated exactly because the 
Coulomb-interaction hamiltonian Hci is quartic in the field operators. In 
the remainder of this article we shall use the approximate* gluon correlator 
1)0(0) evaluated in the ground state ^o^^^ of the harmonic oscillator 
hamiltonian Hho- This introduces no error at high momentum because of 
asymptotic freedom, as is verified in Appendix D where the result obtained 
here is compared with the exact asymptotic expression for u(q) at large q 
calculated using the perturbative renormalization group. It is possible that 
at low momentum the only error is a renormalization of the mass scale. 

In the present section we shall obtain the asymptotic form at high 
momentum of the solution 

* We leave for another occasion to exhibit an improved ground- state wave-function 
obtained by a Bogoliubov transformation. 
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w(q) = u(0)(q) 



(8.1) 



of the integral equation (6.22), in which the kernel is replaced by 
DO(0)i j(k) = [ 5ij - (k2)-lkikj ] DO(0)(k). 

DO(0)k = ( 20)k)-l = 2-1 [ k2 + |x4 w(k) ]-l/2 , (8.2) 

where we have used (7.30) for the frequencies ook of the harmonic 
oscillator hamiltonian Hho- The last equation displays a reciprocal relation 
between the color-Coulomb potential w(k) = u(0)(k), and correlator of 
transverse gluons. The singularity (found in the following section) of w(k) 
at k = suppresses DO(0)(k) at k = 0. 

With this kernel, the integral equation (6.22) for w(q) reads 

w-l(q) = n [2(27i;)3]-l J d3k [ q2 - (q.k)2/k2 ] 

X [ k2 + |i4 w(k) ]-l/2 [ w(k) - w(k+q) ] . (8.3) 

Note that w(q) depends only on the magnitude of its argument. It is 
convenient to write w(q2) instead of w(q), and we have, 

[ q2 w(q2) ]-l = n [2(2jl)2]-l Jo°° dk k2 J.i+1 dz (1 - z2) 

X [ k2 + |i4 w(k2) ]-l/2 [ w(k2) - w(k2+q2+2kqz) ]. (8.4) 
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The Faddeev-Popov operator is positive inside the Gribov horizon, so we 
expect the solution w(k2) of this integral equation to be positive* for all k. 
We shall assume that there exists a unique positive solution. 

We first find the color-Coulomb potential w(q) at momentum large 
compared to the mass scale, q » |i (and also q « a-1). As an Ansatz, 
suppose that at large q, the solution w(q2) differs from the free expression 
w(q2) = l/q2 by a logarithmic power, 

w(q2) = C [ q2 lnP(q2/|i2) (8.5) 
To determine the power p and the coefficient C, we define 

f(q2)^[q2w(q2)]-l, (8.6) 
and apply the operator qd/dq to (8.4). This gives 

q2 f (q2) = . n [2(2ji)2]-l Jo°° dk k2 J.i+1 dz (1 - z2) 

X [ k2 + |i4 w(k2) ]-l/2 (q + kz) w'(k2+q2+2kqz) ]. (8.7) 

We now write k = qx, and take q» [i, and we neglect the term 
|i4w(q2x2) compared to k2 = q2x2. By (8.5) we have, at large q 



* A numerical solution of this equation has been obtained recently which is indeed 
positive [17]. 
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w'(q2) = - C [ (q2)2 lnP(q2/|i2) ]-l [i + p In-l(q2/|i2) ] 

[(q2)2lnP(q2/,i2)]-l, (g.g) 

and 

w'(k2+q2+2kqz) « - C{(q2)2 (x2+l+2xz)2 lnP[q2(x2+l2+2xz)/|l2]}-l 
= - C{(q2)2 (x2+l+2xz)2 lnP(q2/|i2) | - 1 (g.Q) 

This gives asymptotically at large q, 

q2 f (q2) = n C [2(27r)2]-l Jo°° dx x J.i+1 dz (1 - z2) (1 + xz) 

X { (x2+l+2xz)2lnP[q2/|i2]}-l, (g.lQ) 

On the other hand, from (8.6) we have 

q2 f(q2) = C"! p lnP-l(q2/|l2). (8.11) 

Equating the two expressions, we obtain p = 1/2, and 

C-2 = n (27l)-2 jo°° dx X j-i+l dz (1 - z2) (1 + xz) (x2+l+2xz)-2 

C-2 = n (27r)-2 (2/3). (8.12) 
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The color-Coulomb potential at high momentum is now completely 
determined, 

w(q) = (37i;2/n)l/2 [ q2 lnl/2(q/|i) ]-l, q » |i. (8.13) 

Its fourier transform is the color-Coulomb potential at short range in 
position space, 

w(r) = (37t2/n)l/2 (27r)-3 J d3q [ q2 lnl/2(q/^) ]-l exp(iq-x) 

w(r) = (37r2/n)l/2 { 47rrlnl/2[(,xr)-l] r«|i-l. (8.14) 
The color-Coulombic electric field operator 

Eca(x) = - V 1 d3y w(x-y) pa(y). (8.15) 

is expressed in terms of the gradient of this quantity, 

-Viw(r) = (37r2/n)l/2 [ 4n r2 lnl/2(|ir)-l ]-l (xi/r) r « 

(8.16) 

The approximate color-Coulomb interaction energy, 
v(0)(q) = q2w2(q) = q2 [u(0)(q)]2, (8.17) 
is given at high momentum by 
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v(0)(q) = (37r2/n) [ q2 ln(q/|i) q » |io, (8.18) 

and at short range in position space, 

v(0)(r) = (37r2/n) { 4n r ln[ (|xr)-l] r « (8.19) 

It is instructive to write this in the form 

v(0)(r) = gc2(r)/(47rr), (8.20) 
where gc(r) is the Coulomb-gauge running coupling constant 

gc-2(r) = be ln[ (|ir)-l] r « (8.21a) 

be = (37i;2)-l n. (8.21b) 

In Appendix D it is verified that this result agrees with the perturbative 
renormalization group. 

We may use our solution w(k) to evaluate god^-), eq. (6.24). In the 
critical limit, it is given by 

gO-1 = n [3(271)3]-! j d3k [k2 + |i4w(k)]-l/2 w(k). (8.22) 

This integral converges in the infrared for the solution w(k) at low 
momentum which will be obtained in the next section. However, in 
contrast to the integral equation for w(k) which contains the difference 
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[w(k) - w(q)] in the integrand, and is ultraviolet finite, the integral (8.22) 
has a logarithmic ultraviolet divergence. Instead of reverting to lattice 
kinematics, we calculate go(|-i) to leading order in a-1 by inserting a 
momentum cut-off in (8.22) at A = a-1, where a is the lattice spacing. 
This gives 

go(|l)-l = n (27i)-2 (2/3) C lnl/2(A2/,^2)^ (8.23) 
or, by (8.12) 

g0(|^)-2 = be ln(A/|i), A » |i, (8.24) 

where be is the coefficient (8.21) of the running coupling constant. Our 
expressions for unrenormalized coupling constant and for the renormalized 
coupling constant are consistent. 



9. Determination of the Long-Range Color- Coulomb Potential 

We next determine the solution w(q) of the integral equation (8.4) at 
low momentum. As an Ansatz, we suppose that for q « |i, the color- 
Coulomb potential behaves like a power, 

w(q) = B |is-2 q-s. (9.1) 

We change variables in (8.4) according to k = qx, and we now neglect k^ 
compared to [i^ w(k2). This gives 
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B-1 (q/|l)s-2 = Bl/2 (q/|i)3-s/2 n [2(271)2]-! I, 



(9.2) 



where I is independent of q. We thus obtain s = 10/3. For this value of s, I 
is the numerical integral 

I = Jo°° dx xll/3 J.i+1 dz (1 - z2) [ x-10/3 - (x2 + 1 + 2xz)-5/3 ] 
I = 8 r(8/3) r(2/3) r- 1(1 6/3). (9.3) 
This gives 

B-3/2 = n 71-2 r(8/3) r(2/3) r-l(16/3), (9.4) 
u(0)(q) = w(q) = B (|l/q)4/3 q-2, q « ^, (9.5) 

The approximate color-Coulomb potential w(q) is too singular at 
q = to possess a fourier transform, but we may calculate the color- 
Coulomb electric field operator (8.15). If there is an external quark at the 
origin, for example, with color-electric charge Q^, then according to 
(8.15), it produces a static color-electric field 

Ei^ = - Viw(r) Qa = (xi/r) Er (9.6) 

where 

- Viw(r) = - (27t)-3 J d3q iqi exp(iq-x) B |i4/3 q-10/3 (9.7) 



72 



is a finite integral, and the radial color-electric field Er at large r is given 
by 



Er = 3 (167l2)-l 31/2 r(2/3) B |l4/3 r-2/3 j- » ,^-1, (9.8) 
and falls off at large r like r-2/3. 

The color-Coulomb interaction energy at low q is given by 
v(0)(q) = q2Y^2(q)^ or 

v(0)(q) = B2 |i8/3 q-14/3. q « ^. (9.9) 

The strong infrared singularity of this interaction energy provides the 

confining mechanism described in the introduction, eq. (1.18). It 

corresponds to a color interaction potential that grows at large r like 
v(0)(r) ^ r5/3. 

10. Conclusion 

We have developed a Hamiltonian scheme for the minimal lattice 
Coulomb gauge which includes the restriction of the configuration space to 
the fundamental modular region. This restriction implies that Faddeev- 
Popov propagator which transmits the color-Coulomb field is highly 
singular at k = 0. We derived a zero-order hamiltonian Hq, eqs. (1.16) or 
(7.20), that contains an instantaneous color-Coulomb interaction energy 
v(r) between color-charge, and which is the starting point for systematic 
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corrections. An approximate evaluation of v(r) gives an expression v(0)(r) 
that agrees with the perturbative renormahzation-group [r ln(|j,/r)]-l at 
small r, and grows like r5/3 at large r, thus exhibiting the qualitative 
features of a confining theory. 

The question naturally arises to what extent these results are a gauge 
artifact. This term is somewhat of a misnomer since one does indeed have 
the privilege of fixing the gauge. The intended question is to what extent 
the results are an artifact of approximations made after fixing the gauge. 
In particular one may ask whether the rather unexpected r5/3 behavior of 
the approximate color-Coulomb interaction energy v(0)(r) also 
characterizes the exact interaction energy v(r), and if so, and how this is to 
be reconciled with Seller's upper bound on the energy E(r) of an external 
quark pair at separation r, E(r) < Kr, which holds at large r [18]. 

As a general remark it should be noted that the singular behavior in 
the infrared of the Faddeev-Popov propagator, exhibited by eq. (6.9), is an 
exact consequence of the restriction of the configuration space to the 
fundamental modular region. We believe that this is indeed the origin of a 
confining long-range color-Coulomb potential v(r) which appears in the 
zero-order hamiltonian Hq, and that this confining behavior survives the 
corrections due to the perturbation Heff - Hq. 

Whether or not the exact color-Coulomb interaction energy v(r) 
grows at large r like r5/3 as does v(0)(r) requires further investigation to 
determine. The difference v(r) - E(r) is due to the distortion of the gluon 
wave function produced by the external quark pair, and is, in principle, a 
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calculable effect, like vacuum polarization in QED. Moreover in 
Appendix E it is proven that in a world described by Hq, the inequality 
E(r) < v(r) holds. The sign of this inequality and Seller's asymptotic 
bound E(r) < Kr on the energy E(r) of an external quark pair at 
separation r [18], are obviously consistent with E(r) < Kr < r5/3 at 
large r. 

Finally we note that while a long-range instantaneous potential may 
explain confining behavior or an area law for a time-like Wilson loop, it 
says nothing directly about a space-like Wilson loop. In particular, the 
suppression at low momentum of the spatial components of the equal-time 
gluon correlator, (8.2), makes it hard to see in the Coulomb-gauge how 
confining behavior arises for a space-like Wilson loop. This is the same 
difficulty, alluded to in the Introduction, that occurs in the Landau gauge 
for all Wilson loops. That this difficulty is only apparent could be shown 
in the Coulomb gauge by a proof that Lorentz invariance is regained in the 
continuum limit. 
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conversations with Gianfausto Dell'Antonio, Tony Duncan, T. D. Lee, 
Martin Schaden, Dieter Schuette, and Wolfhart Zimmermann. He is 
grateful to the organizers of the Fifth Meeting on Light-Cone Quantization 
and Non-perturbative QCD, Regensburg, Germany, June 20-23, 1995, for 
a stimulating conference that provided the inspiration to start this 
investigation. He is grateful to Professor Lee for sending him a preprint of 
ref. [5] before publication, and to Tony Duncan for bringing ref. [18] to 
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Appendix A: Evaluation of the Kinetic Energy Operator 

The evaluation of T' = 1 - aoK to order ao will be done in a number 
of steps. 

Step 1 . Expansion about saddle-point . 

We wish to evaluate the integral (4.12) for small values of £. The 
variables gx = hx,t that live on vertical links effect a gauge transformation 
on the link variables Ux,t that lives on the time-slice t, for the one- 
plaquette action Sx,i may be written 



As 8 approaches zero, the integrand of (4.12) gets concentrated at the 
minimum of the action of the vertical plaquettes, which occurs at 



namely where B^^ is a gauge-transform of Atr. However a complete gauge- 
fixing has been done on each time- slice, so we conclude that at the 
minimum of the action 



Sx,i = 1 - n-lRe tr [Ut(Ax,itO Ug(Bx,itr) ] . 



(A.l) 



Ug(Btr) = U(Atr), 



(A.2) 



Btr = Atr, 



(A.3) 
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and 

g = 1 (A.4) 

hold. As 8 approaches zero, the integral (4.12) may be evaluated by 

Gaussian expansion around the minimum of the action, in which the 
dominant configurations are of order 

B - A = 0(8l/2) g = 1 + 0(8l/2), (A.5) 

as we shall see. 

We set 

Btr = Atr + 8l/2v, (A.6) 
and expand ^(Bt^) to order 8 

^(Btr) = ^(Atr) + v« 9^(AtO/aAtra 

+ 2-l8 vttvP a2^(Atr)/aAtraaAtrP, (A.7) 
where we have written v instead of v^^ for simplicity. This gives 

(T'^)(Atr) = lo ^(Atr) + a^(Atr)/aAtra 
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+ l2«P a2^(Atr)/aAtraaAtrP, 



(A.8) 



where the Ia(A^^) are coefficients which we now evaluate. 

Step 2. Evaluation of I n. 

Because Iq^ is the only term in this equation that does not contain a 
derivative of ^, consistency with (4.11) requires that, with proper choice 
of normalization constant N, the coefficient Iq be of the form 



where V'(Atr) is a (possible) potential energy contribution to K. By (4.12) 
we have 



where Sx,i is given in (A.l). 

The easiest way to evaluate (A. 10) is to undo the gauge-fixing. By 
(4.5), we have 



lO = N J nx,i dBx,i ll|/(Bx,i)l nx6A[(V-B)x] 

X detM(B) rixdgx exp[- 8 -1 Ix,i Sx,i ] , (A.ll). 



lo(Atr) = [ 1 - ao V'(Atr)], 



(A.9) 



lO = N Ja dBtr 0x4 l¥(Bx,itOI detM(Btr) Oxdgx 



X exp[- 8 -1 Ix,i Sx,i ] , 



(A.IO). 
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where 

Sx,i = 1 - n-lRe tr [Ut(Ax,itr) Ug(Bx,i)] . (A. 12) 

We now use the invariance of the Haar measure dBx,il\|/(Bx,i)l under gauge 
transformation to write 

lo = N 1 nx,i dBx,i l\|/(Bx,i)l exp[- £ -1 Ix,i Sx,i ] 

X rixdgx nx5A[(V-Bg)x]detM(Bg), (A.13) 

where 

Sx,i = 1 - n-lRe tr [Ut(Ax,itr) U(Bx,i)] . (A.14) 

The integral over dgx gives unity by the Faddeev-Popov identity. From 
invariance of the Haar measure under left multiplication we obtain 

lo = N I nx,i dBx,i IV|/(Bx,i)l exp[- £ -1 Lx,i Sx,i ] (A. 15) 

with 

Sx,i = 1 - n-lRe tr U(Bx,i) . 
The last integral is independent of Atr and gives 
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io = i, 



(A. 16) 



which determines the value of the normalization constant N, 



N-1 = nx,i I dBx,i l\|/(Bx,i)l exp{- £ -1 1x4 [ 1 " n'^Re tr U(Bx,i) ] }. 



We conclude that the operator K contains no potential energy term. 

Step 3. Evaluation of vertical plaquette to order £3/2. 

There remains to evaluate the coefficients Ii and I2 defined in (A.8). 
We make the linear change of variable of integration from Rtr to v^^, 
defined by 



and we temporarily write A, B and v instead of Atr, Btr and vtr. We also 
set 



(A. 17) 



and 



V'(Atr) = 0. 



(A. 18) 



Btr = Atr + gl/lytr^ 



(A. 19) 



gx = exp(£l/2cox) 



(A.20) 
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where oOx = t^oOx^ is an (anti-hermitian) element of the Lie algebra. From 
(4.12) we have 

(Il«, I2«,P) = N J dv 11x4 l¥(Bx,i)l detM(B) Ilxdgx 

X exp[- 8 -1 Ix,i Sx,i ] (£l/2v«, 2-l£v«vP) . (A.21), 

where B stands for the right hand side of (A. 19), and the one-plaquette 
action Sx,i is given by (4.13). The limits of integration on v due to the 
restriction to the fundamental modular region are of order 8-1/2, However 
to obtain a Gaussian expansion in powers of £l/2 about the saddle-point at 
V = 0, the limits of the v integration are made infinite. 

Recall that we only need (T'^)(Atr) to order £. Because the last 
factor of (A.21) contains explicit coefficients £l/2 and £, it suffices to 
evaluate the first factor to order £l/2, Xo the order of accuracy required 
we have, 

dgx = do)x. (A.22) 

(We ignore the overall normalization now because it is easier to determine 
it later.) 

In order to evaluate to order £l/2 the first factor of (A.21), which 
contains the expression £-lSx,i, we must evaluate the one-plaquette action 
Sx,i to order £3/2, We write 



81 



Sx,i = 1 - n-1 Re tr expFx,i , 



(A.23) 



where Fx,i is an element of the Lie algebra which, according to (4.13), is 
given by 

expFx,i ^ Ut(Ax,i) gxt U(Bx,i) gy , (A.24) 

where y = x+ei. The leading term in Fx,i is of order £l/2^ and the leading 
term in Sx,i is of order Fxj^, so to evaluate Sx,i to order 8^/2^ it is 
sufficient to evaluate Fx,i to order 8. For this purpose we write 

expFx,i = expF'x,i exp(8l/2o)y) , (A.25) 

where 

expF'x,i = Ut(Ax,i) exp( - 8l/2c0x) U(Bx,i). (A.26) 

We temporarily write A, B and F' for Ax,i, Bx,i and F'x,i. Upon taking 
hermitian conjugates we obtain 

F(B, A, -0)x) = - F(A, B, 0)x). (A.27) 

It is convenient to change argument from A and B to 

Vx,i = 2-1 (Ax,i + Bx,i), (A.28) 



82 



and the variable Vx,i previously introduced, e^^^Vx,! = (Bx,i - Ax,i). We 
have 



expF = Ut(V - 2-l8l/2v) exp( - e^/^Gix) U(V + l-^e^^^w). (A.29) 



so F(V, V, oOx) contains only odd powers of v and oOx- These variables 
appear with coefficient e^^'^, so to obtain F'(V, v, COx) to order 8, it is 
sufficient to expand it up to terms linear in v and oOx, at fixed V. 

For this purpose we write, for an element U(A) of the SU(n) group, 
considered as a function of coordinates Ab, 



where b and c are color indices, the t^ form an anti-hermitian basis of the 
Lie algebra of SU(n), and \|/c(A)dAc is the (Lie-algebra valued) Maurer- 
Cartan one-form of the SU(n) group. The matrix \|/%(A) maps elements 
dAb of the cotangent space of link variables U(A) at A into Lie-algebra 
valued link variables. The Haar measure on the SU(n) group is given by 
l\|/(A)l = det\j/(A). [Explicit expressions for these quantities in the case of 
the SU(2) group are given in eqs. (4.25) to (4.29).] To order e^/^ we have 



and 



F(V, - V, - COx) = - F(V, V, COx), 



(A.30) 



UtdU = \|/c(A) dAc = tb \|/bc(A) dAc, 



(A.31) 



U(V ± 2-l£l/2v) = U(V) [ 1 ± 2-l£l/2xj/d(V)vd ], 



(A.32) 
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so to order 8 we have 

expF' = [ 1 + 2-181/2 x|/d(V) vd ] Ut(V) exp( - e^'^Gix) U(V) 
X [ 1 + 2-181/2 \|/d(V) vd ] 
= 1 + 8l/2 [ \|/d(V)vc - Ut(V)(OxU(V) ], 
which gives, to order 8 

F' = el/2 [ \|/(j(V)vd - Ut(V)C0xU(V) ] . (A.33) 
From (A.25) we obtain, to order 8, 

expF = exp{8l/2 [ \|/d(V)vd - Ut(V)o)xU(V) ] } exp(8l/2o)y) .(A.34) 
By (A.23) we have, for the one-plaquette action, correct to order 8^/2, 

Sx,i = 8 (4n)-l lb { ¥bd(Vx,i)vx,i'^ + (Oyb - [Ut(Vx,i)COxU(Vx,i)]b }2 , 

(A.35) 

where y = x+ei, and we have used tr(tatb) = -2-l5a,b- 

Let D(A) be the matrix that maps Lie-algebra valued site variables co 
into the Lie-algebra valued link variable according to 
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tc [D(A)o)]x,iC ^ 0)y - Ut(Ax,i)o)xU(Ax,i), 



(A.36) 



where y = x+ei. The right hand side is the gauge covariant difference 
between oOy and oOx, where oOx is transported to y by U(Vx,i). The 
operator D(A) is related to the "lattice gauge-covariant derivative D(A), 
introduced in (3.15) by the change of basis 

[D (A)CO]x,ib = V|/ba(Ax,i) [D(A)CO]x/, (A.37) 

where \^^a(^x,0 is the Maurer-Cartan form of the SU(n) group defined in 
(A.31). The adjoint matrices Z)*(A) and \|/*(A) are defined be means of 
the natural inner product of Lie-algebra valued link or site variables, as in 
(3.12). In particular Z)*(A) maps Lie-algebra valued link variables into 
Lie-algebra valued site variables. To prove (A.37), we write 
gx,y = exp(£l/2(j()x y), and we have, to order £l/2 

£l/2 [ COy - Ut(Ax,i)COxU(Ax,i) ] 

= Ut(Ax,i) [ gxtU(Ax,i)gy - U(Ax,i) ] 
= Ut(Ax,i) [ U(Agx,i) - U(Ax,i) ] 
= Ut(Ax,i) {U[ (A + 8l/2D(A)C0)x,i ] - U(Ax,i) }, 
= 8l/2 tb \|/ba(Ax,i) [D(A)0)]x,i^ 
by (3.15) and (A.31). 
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With this result, we may write the one-plaquette action (A. 35), to 
order £3/2^ as 

Sx,i = £ (4n)-l Lb { i|/ba(Vx,i) Vx,i^ + [£>(Vx,i)0)]x,i^ , (A.38a) 
or, in matrix notation 

Sx,i = £ (4n)-l ( \|/(V) V + D(V)co, \|/(V) v + Z)(V)co ) , (A.SSa) 

where V = 2-l(A+B), and v is a vector with independent components y*^, 
that represents a generic configuration defined on horizontal links of the 
lattice. This expression has a simple geometric interpretation. Recall that 

ds2 = lb [ ^K(A) dAc ]2 . (A.39) 

is the line element at A of the SU(n) group that is invariant under left and 
right multiplication. With B-A = e^^^v, and g = exp(e^^^(d), we see that 
to order £3/2, the one-plaquette action is the line element of the difference 
of the horizontal links B-A, plus the gauge covariant derivative of the 
gauge-transformation g defined by the vertical links, evaluated at the mean 
V = 2-l(B+A). This expression is a lattice analog of the continuum 
expression Foi^ = QoM - Di(A)Ao)^. 

Step 4. Integration over vertical links, and lattice form of Gauss's Law 
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We now integrate out the variables oOx associated with the vertical 
links by Gaussian integration, with the result 

J dco exp[- 8 -llx,i Sx,i] 

= det-l/2[/)*(V)Z)(V)] exp[ - (4n)-l (v, G(V)v) ]. (A.40) 

The matrix G(V) defines a metric tensor on transverse link configurations 
which is the lattice analog of the gauge-invariant metric introduced in 
continuum theory by Babelon and Viallet [19]. Its explicit expression is 

G(V) = Ptr \|/*(V) [1-D (D*DyW* ] \|/(V) Ptr, (A.41) 

where D = D(V) and Z)* = Z)*(V). Here Pt^ is the projector onto trans- 
verse configurations 

Ptr = 1 - V(V*V)-lV*. (A.42) 

It appears on the left and right of G, because in (A.41), v is transverse, 
V *v = 0, so the metric tensor G(V) in (A.40) is defined only on 
transverse configurations. Each factor in (A.41) is a matrix that acts on 
link variables. The middle factor in G(V), 

PG = 1 - D (Z)*Z))-lZ)* (A.43) 
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is the projector onto the subspace that satisfies Gauss's law, D*(V)PG = 0, 
and in consequence, the electric field operator defined below will satisfy 
the lattice form of Gauss's law exactly. 

Step 5. Integration over horizontal links 

With this result, we obtain from (A.21) to order 8, 

(II", l2a,P) = N 1 dv det\|/(B) detM(B) det-l/2[Z)*(V)Z)(V)] 

X exp[ - (4n)-l (v, G(V)v) ] (£l/2va, 2-l8v«vP) , (A.44), 

where B = A + el/^v, and V = A + 2-^£^/^\, det\|/(A) = nx,idet\|/(Ax,i). 
Because I2 contains an explicit factor of 8, we may set 8 to zero 
everywhere else. This gives 

l2a,P = Iq(0) n8 [G-1(A)]«'P, (A.45) 

where 

Io(0) = C det\|/(A) det[D*(A)V] det-l/2[Z)*(A)D(A)] det-l/2G(A), 

(A.46) 

M(A) = D*(A)V, and C is a constant which we shall determine shortly. 
Here Iq^^^ is the value of Iq, to zeroth order in 8, which one obtains if Iq is 
evaluated by direct integration, as we have done for I2. However by 
(A. 16) we have, Io = l, which is exact for all 8 and A, so also 
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io(0) = i. 



(A.47) 



This determines the value of the constant C, which is found by setting 

A = 0, namely C = det-l/2(V* V). It also yields the value of detG(A), 
namely 

detl/2G(A) = det\|/(A) det(D*(A)V) 

X det-l/2[Z)*(A)Z)(A)] det-l/2(V*V). (A.48) 

We thus obtain for 

I2«'P = ne [G-1(A)]«'P. (A.49) 

There remains to evaluate Ii^^. Expression (A.44) for li^ contains 
an explicit factor of £l/2 in the right hand factor, so it is sufficient to 
evaluate the first factor to order 8 1/2, Therefore in the first factor, where 
B = A+£l/2v, and V = A+l-^e^^^w, we may set 

det\|/(B) detM(B) = ( 1 + £1/2 yP d/dA^ ) [det\|/(A) detM(A)] 

(A.50a) 

and 

det-l/2[D*(V)Z)(V)] exp[ - (4n)-l (v, G(V)v) ] 
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= (1 + 2-181/2 vP a/aAP) {det-l/2[Z)*(A)Z)(A)] 

X exp[ -(4n)-l (v, G(A)v)]}. (A.50b) 

This gives, with 

o(A) = det\|/(A) detM(A) (A.5 1) 

II" = 2n8 a-l(A) [G-l(A)]a.P do(A)/dA^ 

+ ne a(A) a/aAP{ a-l(A) [G-l(A)]a,P } 

II" = n£ o-l(A) a/aAP {o(A) [G-1(A)]«'P}, (A.52) 

and we have used (A.48) 

Upon substituting these results into (A. 8), we obtain finally, with 
8 = go^ ao (2na)-l, 

(T'^)(A) = ^(Atr) - ao K^(A), (A.53) 
where K is the kinetic energy operator 

K = - go2 (2a)-l a-l(A) d/dA^ a(A) [G-l(A)]a,P d/dA^. (A.54) 
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Step 6. Lattice color-electric field operator and lattice color-Coulomb 
operator 

The explicit form of G"1(A) is given by 

G-1 = [ 1 - D(V*D)-lV*] J* J [ 1 - V(D*V)-lD*], (A.55) 

where D = D(A), and V*D = D*V=M, where M = M(A) is the 
Faddeev-Popov operator. Here J = J(A) is a matrix which is diagonal on 
links that is inverse to the Maurer-Cartan form, namely on each link it 
satisfies 

JaKAx,i) ¥^c(Ax,i) = 6cb (A.56) 
The is^(Ax,0 satisfy the Lie differential equations of the SU(n) group, 

[Ja, Jb] = fabcJc, (A.57) 

where 

Ja ^ Jab(Ax,i)a/aAx,ib. (A.58) 
To prove (A.55), we shall verify that 

G Q = Ptr, (A.59) 
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where Ptr is the unit operator on the space of transverse link variables. 
Here Q is the expression on the right hand side of (A.55), and G is given in 
(A.41). We obviously have V*Q = QV = 0, so 

Ptr Q = Q, (A.60) 

and moreover 

\|/ Ptr Q = \|/ Q = [ 1 - Z)(V*D)-1 V*J*] J [ 1 - V(D*V)-lD*], 
where D = \|/D. This gives, as required, 

G Q = ptr \|/* [ 1 - Z) (Z)*D)-lZ)* ] [ 1 - Z)(V*D)-1 V*J* ] 

X J [ 1 - V(D*V)-lD*] 
= ptr xj/* [ 1 - D ] J [ 1 - V(D*V)-lD*] 

= ptr \|/* [ J - Z) (Z)*Z))-lD* ] [ 1 - V(D*V)-lD*] 
= ptr \|/* J [ 1 - V(D*V)-lD*] = ptr [ 1 - V(D*V)-lD*] = Ptr, 
where we have used Z)*J = D*. 

Appendix B: Explicit form of Transformed Hamiltonian 
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After transformation (5.15), the hamiltonian is given by the 
quadratic form 

(<!>', Heff' = (<!>•, K' + (<!>', (Vm + YoG) ^O, (B.l) 

where Vm is the magnetic field energy of (4.16), G is given in (5.1), and 
the transformed kinetic energy operator is given by 

(VF-, K'^)') = 1 dAtr [go2(2a)-l] Ix,i (£x,i,a^' - (2a)-l [£x,i,a, CT]0')* 
X (£x,i,a^' - (2o)-1 [ £x,i,a, CJ] W), (B.2a) 

where the operator 

Ex,a,i = £x,a,i«( Atr) i 3/3 Atr,a (B .2b) 

is given in (4.17) and (4.18). With 

Ex,a,i^ ^ i a/aAtr,a£x,a,i«(AtO, (B.2c) 

the transformed hamiltonian operator acts according to 

Heff' = K'^' + (Vm + YoG) (B.3) 

where 
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K'^' = [go2(2a)-l] Ix,i {E^,!^"^ £x,i^ - [^x,i,a« (2a)-iaa/aAtr,a ]2 
+ 3 [£x,i,a« £x,i,aP (2a)-iaa/aAtr>P] /aAtr>a } vp'. (B.4) 
The last two terms represent an effective potential energy. 



Appendix C. Properties of the Faddeev-Popov Propagator 

We write the Faddeev-Popov matrix as 

M = Mo + Mi, (C.l) 

where 

Mo = V*V (C.2) 

Ml = go B*V = go V*B (C.3) 
and the operator B is defined by 

go B ^ D - V, (C.4) 
where D is the lattice gauge-covariant derivative of sect. 3. We have 

M-1 = Mo-1 - Mo-1 Ml Mo-1 + Mo-1 MiM-IMiMq"!. (C.5) 
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By the definition of D 

(Diaccoc)x = Gx,i^^ViCOC + go fabcAx^baiCOx^, (C.6) 

where aiCOx = 2-l[(0(x+ei) + 0)(x)], the matrix M"! has the particular 
structure 

M-1 = Mq-I + Mq-I V*FVMo-1 - go Mo-1 V*A Mq-^. (C.7) 

Here the action of the operator A is defined by the last term of (C.6), and F 
is defined by 

Fx,i,a;y,j,b ^ gO^ B(x)ia,c B(y)jb,d (M-l)x,c;y,d 

- (Gx,i^'b - 6a,b) 5x,y 5i,j. (C.8) 

The explicit expression for Gx,i^'b is given in (3.7). 

Consider the related matrix E defined by 

Ex,i,a;y,j,b = Eix,i,a;y,j,b + E2x,i,a;y,j,b (C.9) 

where 

Elx,i,a;y,j,b ^ D(x)ia,c D(y)jb,d (M-l)x,c;y,d. (C.IO) 

E2x,i,a;y,j,b = - Gx^^'b 5x,y 5i,j. (C.ll) 
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Because D = goB + V, the matrices E and F are related by 

Fx,i,a;yJ,b = [ Ex,i,a;y,j,b + 6x,y 6i,j 6a,b ]TT (C.12) 

where the superscript TT means the part that is transverse on both indices. 
The matrix E is related to the horizon function G by 

G = Lx,y,i,a Ex,i,a;y,i,a- (C.13) 

The propagator C, defined in (6.1) may be written 

C = Mo-l +M0-I V*<F> VMq-I, (C.14) 

where we have used <A) = 0, which follows from global color or lattice 
symmetry. In terms of its fourier transform (6.2), this reads 

Ce = + Ae-^ qe,iKe,i,j qe,j (C.i5) 



where 



Ke,i,j 6a,b ^ v-1 Ix,y exp[-ie-(x + 2-lei - y - 2-lej)] <Fx,i,a;y,j,b>. 

(C.16) 



We also define 
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Le,i,j 6a,b ^ v-1 Ix,y exp[-ie-(x + 2-lei - y - 2-lej)] <Ex,i,a;y,j,b>, 

(C.17) 

which has the decomposition 

Le,ij = Lie,ij + L2e,ij- (C.i8) 

corresponding to the decomposition E = Ei + E2. By (C.12), Ke^j and 
Le4j are related by 

Ke,ijTT= [ Le,ij + 5ij ]TT (C.19) 

We now estabhsh a number of properties of Le^j. By (C.l 1) and by 
lattice and global color symmetry we have 

L2e,iJ = - b 6ij, (C.20) 

where 

b 5a,c = <Gx,ia,c> = 5a,c <n-l Re tr Ux,i>, (C.21) 
by (3.7). The matrix Ei satisfies 

V*Eix,a;yj,b = D(y)jb,d 5x,y 6a,d, (C.22) 
and thus 
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Li qe,i Lie,iJ = b qe,j, 



(C.23) 



where we have again used <A) = 0. This gives 

Lie,iJ = b 5ij + Se,ijTT, (C.24) 
where SejjTT transverse on both indices. We have 

Le,ij = Se,i,jTT, (C.25) 

so L04J is transverse on both indices. From (C.13), the horizon condition, 
<G) = 0, is equivalent to the statement 

Li Le,i,iie=o = 0. (C.26) 



So far the volume of the lattice is finite. We now take the infinite 
volume limit, so 0i becomes a continuous angular variable, and we write 
Lij(0), and similarly for other 0-dependent quantities. We assume that 
global color symmetry and cubic symmetry of the lattice are preserved in 
the infinite-volume limit. We also assume that the leading terms at 0i = 
have the invariance of the full (3 -dimensional) Euclidean group and not 
just the symmetries of the cubic lattice group. This means that Lij(0) and 
Kij(0) are Euclidean-invariant functions of q(0) at small 0. Then (C.25) 
implies that in the neighborhood of 0i = 0, the tensor Lij(q) is of the form 

Lij(q) - a(q2) ( 6ij - qi qj /q2 ), (C.27) 
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because this is the only transverse tensor with the stated symmetries. By 
eq. (C.26), the horizon condition is the statement 



a(0) = 0, (C.28) 

and we have Lij(0)= 0. According to (C.19), this is equivalent to the 
statement 

[ K(0)ij - 6ij ]TT = 0, (C.29) 
so the horizon condition has been expressed as a property of C(0). 



We wish to express C(0) in terms of its self-energy. Let the matrix 
E(9) with elements I^i,j(G) be defined by 

zee) = K(0) [ 1 + P(e) K(e) (c.30) 

where P(6) is the projector onto the longitudinal subspace, 

Pij(e)^qi(0)qj(e)/q2(e). (C.31) 

It has the expansion 

ZijCe) = Kij(e) - Ki,k(e) qkO) [q2(e)]-l qiO) KijO) + ... . 

(C.32) 
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Equation (C.30) may be inverted to give 



K(0) = 2:(e) [ 1 - P(e) 2:(e) 



(C.33) 



If this is substituted into eq. (C.15) for C(0) one obtains (6.3). 

The expansion (C.32) shows that the transverse parts of Kij(9) and 
Zij(0) are equal at 6 = 0, 



and we conclude that the horizon condition in the form (C.29) is equivalent 
to (6.6). 

In expression (C.32) for I^e,i,j, the leading singularity of Kij(0) at 
= 0, is canceled by the remaining terms. For Kij(0), we use (C.16) and 
(C.8), and obtain from (C.32) 



[Eij(0)]TT = [Kij(0)]TT , 



(C.34) 



^QXi = [(n2-l)V]-l Ix,y exp[-i0-(x + 2-lei - y - 2-lej)] 



X go2 <B(x)ia,c B(y)ja,d (M-l)x,c;y,d> + ... . (C.35) 



This gives 



^ [(n2-l)V3]-l Ix,y,e',e" exp[-i0-(x + 2-lei - y - 2-lej)] 
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X go2 <(M-l)e',c;e",d B(x)ia,c B(y)ja,d > exp[i(e'-x - e"-y) ] + ..., 

(C.36) 

where (M-l)e',c;0",d is defined in (6.11). We use (C.4) and (C.6) for 
B(x)-a,c^ and write explicitly only the leading term in go, to obtain (6.10). 

Appendix D. Comparison with the Renormalization-group and 
Perturbative Calculations 

The energy E(r) of a pair of external quarks at separation r is a 
gauge-invariant quantity of the form 

E(r) = (47rr)-l grg2(r) (Qi^ Q2a), (D.l) 

where grg(r) is the renormalized running coupling constant of QCD. At 
small separation r it is given by [20], 

grg-2(r) = brg ln[ (rAqcD)"! ] , r « AqcD"! (D.2) 

brg = (247i;2)-l (1 In - 2nf ) , (D.3) 

where nf is the number of quark flavors, and Aqcd is a hadronic mass. 
Expressions (8.21) and (D.2) exhibit the same analytic dependence on r, 
but the numerical coefficients be and brg are different. 
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To understand this difference, it is sufficient to consider the 
perturbative calculation of E(r) up to one-loop order, because the 
renormalization-group coefficient brg is determined by the one-loop 
perturbative calculation. The perturbative calculation to one-loop in the 
Coulomb gauge is presented in [20], with the result 

E(r) = (47ir)-l gp2(r) <Qia Q2^), (DA) 

where 

gp2(r) = go2 + go4 brg In(rA) + ... (D.5) 

and A is an ultraviolet cut-off. The go^ term is the classical Coulomb 
energy, and the go"^ term is the one-loop contribution, and brg is the 
renormalization-group coefficient. In the Coulomb gauge, it results from 
the sum of two contributions [20], 

brg = 5 + 6' (D.6) 

where 

6 = (2ji2)-l n (D.7) 

6' = - (247t;2)-1 ( n + 2nf ), (D.8) 

whose physical origin we now explain. Both terms result from the part of 
the Coulomb hamiltonian given by 
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Hi = 2-1 L 



Px^ [ go (M-l)ab( - V2 )-l go (M-l)bc]x,y pyC, 



(D.9) 



where p = pdyn + Pext- The contribution 5' results from the second-order 
quantum-mechanical perturbative formula 



where to order go^, one makes the replacements M-1(-V2)M-1 — > (-V2)-l. 
This contribution is purely negative and corresponds to the familiar 
vacuum polarization or screening of QED. Thus 6' should not contribute 
to the static color-Coulomb coefficient be, because it comes from a 
distortion of the wave-function caused by the presence of external charges, 
and this is not part of what we have called the static color-Coulomb 
interaction. 

The contribution 5, which is positive and represents anti-screening, 
comes from treating Hi in first-order quantum-mechanical perturbation 
theory 



Ed) = 2-1 Ix,y (01 Pext,x^ [(goM-l)ab( - V2 ) (goM-l)bc]x^y p^^^ yC |0>. 



E(2) = -Inl<0l Hiln) |2 (En - Eo)-l, 



(D.IO) 



E(l) = <0I HilO), 



(D.ll) 



evaluated to order go^. We have 
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Ed) = <0I [ (goM-l)ab( .V2) (goM-l)bc ], q I0> (Qi^ Q20). (D.12) 



The Faddeev-Popov operator, M = - - ggA,, where X = A V, has the 
expansion in powers of go 

M-l( -V2) M-1 = ( - V2)-l + 2 ( - V2)-l go^l ( - V2)-l 

+ 3 ( - V2)-l go^i ( - V2)-l go^i ( - V2)-l + ... (D.13) 

The go^ contribution comes from the last term, with coefficient 3, which 
gives the contribution 5 to brg. 

We may evaluate separately the contribution to E(l) that comes from 
the vacuum intermediate state, and the part that comes from excited 
intermediate states, 

<0I M-l( -V2) M-1 10) = <0I M-1 10) ( -V2) <0I M-1 10) 

+ Ln^O <0I M-1 ln>( -V2) <nl M-1 I0>. (D.14) 

The two parts contribute additively to 6 and brg, 

bi.g = 5 + 5' = 5v-^6e-H5', (D.15) 

Only the first term in (D.14) corresponds to what we have called the color- 
Coulomb interaction energy, because it is obtained by the mean-field 
substitution M-1 <M-1), so consistency requires be = 5v. 
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The perturbative expansion of the vacuum-intermediate- state 
contribution to (D.14) is given by 

<M-1>(-V2)<M-1> = (- V2)-l 

+ 2 < ( - V2)-l go^l ( - V2)-l go^l ( - V2)-l > + ... (D.16) 

The last term has coefficient 2 as compared to the coefficient 3 in (D.13), 
and so, by (D.7), 

6v = (2/3) 5 = (37r2)-l n, (D.17) 

and indeed, by comparison with (8.21) one has, as required 

bc = 6v. (D.18) 

We conclude that the lowest order calculation in the present 
approach gives directly the renormalization-group result for the static 
color-Coulombic interaction energy. A renormalized running coupling 
constant gc(r) has been obtained in the present approach without any 
renormalization having been done. The mass |l corresponds to Aqcd- 

Qualitatively, be is the largest contribution to brg, being 8/1 1 of the 
gluonic contribution to brg. For nf = 6 flavors of quarks and for n = 3 
of SU(3), the ratio is bc/brg = 8/7. Although the corrections to be are 
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numerically significant (and calculable), they do not produce a qualitative 
change. 

Appendix E: Lower bound on quark-anti-quark potential 

Consider the hamiltonian Ho,tot, eq. (7.34), which describes an 
external quark and antiquark at lattice sites and r in interaction with 
dynamical gluons. (Dynamical quarks would screen the interaction of 
interest.) The color-charge density operator of the external quarks is given 



where t(i)^ and t(i)^* are anti-hermitian generators of SU(n) in the 
fundamental representation and its complex conjugate representation 
respectively. The hamiltonian (7.34) has the form 



by 



Pext,x^ = i ( t(l)a 5x,0 + t(2)a* 6x,r ), 



(E.l) 



Ho,tot — Hgiue + Hcross-term + Hext? 



(E.2a) 



where Hgiue = Hq, eq. (7.20), 



Hext = [2 (271)3]- 1 J d3q Pext^(q) v(q) Pext^q) , 



(E.2b) 



and H, 



cross-terai contains the terms pext v pgiue- 



As trial wave-function we choose 
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(E.3) 



where t(l) and t(2) act on the first and second indices respectively, and 
is the vacuum or ground state of Hgiue, so Hgiue^O = 0. The quantum- 
mechanical energy E(r) of the ground state of the system with the external 
quarks satisfies the bound E(r) < F(r), where 

F(r) = C¥, Ho,tot (E.4) 

For the wave-function (E.3), Hcross-term gives vanishing contribution to 
F(r), and we obtain 

F(r)^(^,Hext^). (E.5) 

The color-charge density operator (E.l) acts on the wave-function (E.3) 
according to 

(Pext,x^ ^)a,p = i n-1/2 (t)ap ( 5^ q . 5^ ) vj/q, (E.6) 
because [(ta)a,p]* = - (ta)p,a, which reads, in momentum space, 

(pextHq) ^)a,p = i n-1/2 (t)ap [ 1 . exp(iq r) ] ^0 • (E.7) 
We obtain for F(r), 

F(r) = [2(271)3]-! n-1 (n2-l) j d3q v(q) [ 1 - cos(q-r) ], (E.8) 
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where we have used tr(-t^ t^) = 2-l(n2-l). For large Irl, we may replace 
v(q) by its asymptotic form at low momentum v^s^q), xhis gives for the 
energy of a quark-antiquark pair at large separation r, 

Fas(r) = [2(271)3]-! n"! (n2-l) j d3q vas(q) [ 1 - cos(q-r) ]. (E.9) 

For the asymptotic form which is a power law, vas(q) = c/qP we 

obtain 

Fas(r) = (27i;)-2 n"! (n2-l) C rP-3Jo°° dx x2-p [ 1 - x'lsinx ], 
which converges for 3 < p < 5, with the result 

Fas(r) = (27r)-2 n"! (n2-l) C rP-3 

X [(p-2)(p-3)(p-4)]-l sin(7ip/2) r(5-p) . (E.IO) 

If vas(q) has the form (9.9), with p = 14/3, one has 

Fas(r)~r5/3. (E.ll) 

This expression for F^s and the upper bound E^^(t) < F^^(t) are consistent 
with a linearly rising quark- anti-quark energy E^s^i-) = Kr. 
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Figure Captions 

Fig. 1. Plots of f(cx) = N-l(<G)a + oc-1) against oc, see eqs. (2.22), (2.24). 

(i) The number of degrees of freedom N is finite. 

(ii) As N ^ oo, ttcr 0. This is case (a), the perturbative phase. 

(iii) As N ^ oo, CXcr remains finite. This is case (b), the horizon phase. 
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